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ABSTRACT 

This  work  is  an  attempt  at  determining  the  optimum 
neutral  axis  location  in  the  midship  section  of  a  vessel, 
using  as  a  model  for  the  calculations  a  box  girder,  subject 
to  several  different  conditions  of  sagging  and  hogging 
bending  moments-  several  values  of  design  stress  are  also 
considered  and  in  the  end  what  happens  when  lateral  hydrosta- 
tic load  is  considered  along  with  different  values  assigned 
to  the  design  stress  in  tension  and  in  compression  is  also 
analysed. 
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LIST  OF  SYMBOLS 

a  -  distance  between  transverse  stiffening  members 

AL  -  length  of  box  girder  or  vessel 

B  -  beam 

BOT  -  thickness  of  the  plating  on  the  bottom 

BN  -  bending  moment 

D  -  depth 

DIBOT  -  distance  of  the  bottom  from  the  neutral  axis 

DITOP  -  distance  of  the  top  from  the  neutral  axis 

DSC  -  design  stress  in  compression 

DST  -  design  stress  in  tension;  also  design  stress  without 

discrimination  when  DSC=DST 

H  -  head  of  sea  water  (  lateral  load  ) 

HB  -  head  of  sea  water  on  the  bottom 

HM  -  hogging  bending  moment 

HT  -  head  of  sea  water  on  the  top 

I  -  moment  of  inertia  of  the  cross  section  with  respect 

to  the  neutral  axis  (  GMI  in  the  program  ) 

LBS  -  longitudinal  bending  stress 

S  -  spacing  of  longitudinal  stiffening  members 

SBCRIT  -  critical  stress  on  bottom 

SBH  -  stress  on  bottom  due  to  hogging 

SBS  -  stress  on  bottom  due  to  sagging 

SM  -  sagging  bending  moment 

STCRIT  -  critical  stress  on  top 

STH  -  stress  on  top  due  to  hogging 

STR  -  stress 


STR1  -  primary  stress  =  LBS 

STH2  -  secondary  stress 

STR3  -  tertiary  stress 

STR3B  -  tertiary  stress  on  bottom 

STR3T  -  tertiary  stress  on  top 

STRCR  -  critical  stress 

STS  -  stress  on  top  due  to  sagging 

t  -  thickness  of  theplate 

T  -  draft 

TOP  -  thickness  of  the  plating  on  the  top 

TS  -  thickness  of  the  side  plating 

TSTH  -  total  stress  on  top  in  hogging 

TSTS  -  total  stress  on  top  in  sagging 

TSBH  -  total  stress  on  bottom  in  hogging 

TSBS  -  total  stress  on  bottom  in  sagging 


INTRODUCTION 

It  is  frequently  found   in  the  literature  related  with 
the  strenght  of  ships  such  statements  as: 

"In  general  the  neutral  axis  is  expected  to  be  at  40$  of 
the  depth  (D)  above  the  base  line  ...  a  preliminary  estimate  of 
the  thickness  required  for  the  bottom  plating  (BOT)  can  be 
obtained  from 


0 . 4D        3  / 

1.875( )2240  \/  AL  =■ 

0.6D        V 


40300 


1.75 
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■(■ 
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950    BOT 


where  the  Montgomerie1 s  expression  is  being  considered  to 
determine  the  critical  stress  (SBCRIT  psi)  and  the  short  edges 
are  assumed  far  apart  and  simply  supported. 


SBCRIT=- 


40300 


1     S 

1+ ( ) 

950   BOT 


1.75 


S>4a 

a-  distance  between  transverse  stiffening  members  (inches) 
1.875  accounts  for  a  factor  of  safety 

Tobin's  expression  was  considered  in  the  above  formula 
to  give  the  maximum  permissible  bending  stress  (LBS,  psi) 


LBS=  2240 


AL 
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We  see  thus  that  this  procedure  allows  for  an  estimate  of 
the  bottom  thickness  to  be  obtained  after  knowing  the  length  and 
assuming  a  longitudinal  frame  spacing.  What  may  be  questioned 
here  is  the  reason  why  one  assumes  the  neutral  axis  0.4D  above 
the  base  line  and  also  the  reasons  behind  the  choice  of  a  frame 
spacing. 

In  the  following  a  procedure  to  estimate  the  scantlings 
of  longitudinal  members  is  described;  this  procedure  can  also 
be  found  frequently  in  the  literature  dealing  with  the  subject. 

. . . "For  a  particular  set  of  principal  dimensions  ,  general 
arrangement  and  frame  spacing,  plate  panel  thicknesses  for  one 
after  another  of  the  amidships  longitudinally  continuous  elements 
are  selected  for  an  assumed  position  of  the  neutral  axis  at  kQ% 
of  the  depth  above  the  keel.  The  composite  midship  section  assem- 
blage is  thentested  for  the  adequacy  of  its  section  modulus  when 
ranged  alongside  some  standard  value,  such  as  given  by  the 
Load  Line  Regulations  or  some  function  of  it.  Thus  the  modulus 
is  judged  sufficient  to  limit  ship  bending  stresses  to  the  value 
previously  set  and  incorporated  into  the  design  equations. 

For  all  but  the  deepest  ships  the  modulus  will  undoubtedly 
need  to  be  increased  and  a  small  arbitrary  increment  is  added 
to  the  upper  deck  or  bottom  whichever  happens  still  to  be  farther 
from  the  true  neutral  axis  location  found  after  the  first  design 
cycle.  The  shift  of  neutral  axis  resulting  from  the  scantling 
increment  requires  a  compensating  increment  in  the  opposite  flan- 


ge  with  the  obvious  outcome  that  the  neutral  axis  is  shifted  in 
reverse  with  attendant  effect  on  the  scantling  chosen  first;  a 
cyclic  process  is  then  begun.  More  trials  will  be  necessary 
before  the  match  of  available  to  required  section  modulus  is 
achieved. " 

From  this  transcript,  we  may  again  note   the  initial 
statement   "...for  a  particular  set  of  principal  dimensions, 
general  arrangement  and  frame  spacing. . . "  and  question  whether 
or  not  the  frame  spacing  should  be  chosen  as  an  input  or  rather 
be  left  as  a  variable  and  be  obtained  as  a  result  of  a  structural 
design  procedure.  It  is  here  aknowledged  though  that  in  these 
procedures  an  acceptable  solution  can  be  worked  out  and  each 
new  cycle  is  an  improvement  over  the  precedent  ones.  What  may 
be  questionable  though  is  whether  the  first  cycle,  due  to  the 
assumptions  made  at  its  beginning,  has  or  not  a  strong  influence 
in  the  outcome  leading  to  a  solution  among  the  various  possible 
ones  which  in  the  end  does  not  lead  to"the  optimum  design" 
whatever  its  definition  may  be. 

It  seems  that  this  trial  procedure  can  be  looked  at  as  a 
"weakness"  of  the  structural  design  procedure  for  the  determina- 
tion of  scantlings  and  it  is  perhaps  reasonable  to  make  some  effort 
in  order  to  work  out  a  process  that  could  give  the  scantlings 
in  a  direct  manner  and  be  at  the  same  time  in  agreement  with 
what  can  be  considered  "the  optimum  design",  (usually  based  on 
economic  considerations). 
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Nowadays  much  concern  is  given  to  economics  in  ship  design 
and  naturally  the  strength  calculations  of  a  ship  are  also  made 
with  economical  optimization  in  mind.  With  this  respect,  we  can 
realize  that  in  ships  structures  there  is  a  desirability  to  mi- 
nimize structural  weight  and  some  define  structural  efficiency 
as  being  the  least  weight  structure  capable  of  resisting  the 
loads  for  which  they  are  designed. 

It  should  be  clear  that  the  cost  of  construction  and  of 
maintenance  and  repair  are  inter  related  factors  together  with 
least  weight,  to  be  considered  when  bringing  economics  into 
ship  structural  design.  Also  other  matters  such  as  accessibility, 
weldability  and  availability  of  material  must  be  considered  as 
well. 

In  the  following  sections  an  attempt  is  made  to  obtain 
without  any  interactions  the  scantlings  for  the  plating  of  a 
box  girder.  It  is  hoped  that  through  this  simple  model  some 
helpful  conclusions  can  be  made  in  view  of  their  application 
to  a  vessel's  structure.  Each  of  the  different  situations  ana- 
lized  will  be  exhaustively  explained  in  the  following  sections. 
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PRELIMINARY  CALCULATIONS 


The  model  used  in  this  section  and  also  throughout  this 
work  with  only  very  small  additional  considerations  is  a  box 
girder  which  cross  section  is  pictured  in  figure  1. 


Fig.  1 


As  will  be  seen  in  the  calculations  of  the  following 
sections,  the  dimensional  parameters  used  are  of  the  order  of 
magnitude  of  those  considered  in  ship  design. 

This  model  was  chosen  due  to  its  simplicity  for  the 
performance  of  the  following  calculations  as  can  be  understood 
immediately  from  its  geometrically  symmetrical  cross  section. 

The  objective  of  this  section  is  to  determine  the  neces- 
sary formulation  to  be  used  in  future  sections. 

The  formulation  dealt  with  in  this  section  is  related 
only  to  the  particular  model  chosen  and  the  geometry  of  its 
cross  section. 
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Z(A*d) 

DIBOT= *12 

Z  A 

2ft*d) 

DIBOT'= 

HA 


2  A=B*T0P+2*TS*D+B*B0T=2*D*TS+B*( TOP+BOT) 

DIBOT  =Neutral  axis  location  above  bottom  (in) 
DIBOT'=    "      "      "       "      "     (ft) 

BOT  =Moment  of  inertia  with  respect  to  the  bottom 

I  =   "     "    "     "      "    "    "   neutral  axis 

Z(A*d)=  ( D/l  2  )  *  (  B*TOP+D*TS  ) 

I=Z(A*d*d)f  2  I- (  Z(A*d» 2/ (  2  A ) 

DI BOT=D* ( B*TOP+D*TS ) / ( 2*D*TS+B* ( TOP+BOT ) ) 
This  formula  will  give  DIBOT  in  inches  if  D  and  B  are 
expressed  in  inches;  it  will  give  DIBOT  in  feet  in  the  event 
that  D  and  B  are  expressed  in  feet. 

DITOP=D-DIBOT 

DITOP=D* ( B*BOT+U*TS ) X2*D*TS+B ( TOP+BOT ) ) 

DITOP  in  units  of  inches  or  feet  as  explained  for  DIBOT 
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The  following  formula  will  give  I  in  units  of  (in*ft)2 
for  B  and  D  given  in  inches. 

( D/l 2 ) 2* ( 3*B*B*T0P*B0T+ ( D*TS ) 2+2*B*D*TS* ( TOP+BOT ) ) 


1  = 


3* ( 2*D*TS+B* ( TOP+BOT ) ) 


2 
Next  formula  gives  I  in  units  of  (in*ft)   for  B  and  D 


given  in  feet 


I=- 


2,  2 

4*D*D*  (  3*B*T0P*B0T+  ( D*TS )  +2*B*D*TS*  ( TOP+BOT  )  ) 


2*D*TS+B*( TOP+BOT) 

These  expressions  will  be  used  in  most  of  the  following 
sections  with  B  and  D  expressed  in  feet. 

NOTE:  TS  which  is  defined  as  the  thickness  of  the 
side  plating  is  considered  throughout  this  work  to  be  uniform 
and  will  always  be  equal  to  the  larger  of  the  thicknesses  obtained 
for  the  top  and  for  the  bottom  plating. 
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neutral 
surface 


FORMULAS  USED  IN  THE  ANALYSIS 

BEAM  EQUATION 

Assuming  pure  bending,  the  stresses  in  the  cross  section 
of  the  box  girder  are  given  by 

BM*y  r 

STR= 

I  " 

Fig.  2 

where  y  is  the  distance  from  the  neutral  axis  to  the  point 

where  the  stress  is  being  determined. 

The  neutral  surface  indicated  in  fig.  2  is  composed  of  all 

the  points  which  are  not  subjected  to  any  stress  in  pure  bending. 

BRYAN  EQUATION  -  Plate  under  inplane  load 


(fl)2*E*(t)2 


STRCR=K*' 


12*(1-(P0I)2)*(S) 


a 

fc . -. 

*■  S 

»-  4 

fc-l — » 


psi . 


When  a/S>2  then  K=4 

STRCR  given  in  psi  for  t  andS  expressed  in  inches  and  E  in 

#=3.1^16 
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PLATE  UNDER  LATERAL  LOAD 


S   ?    1 


STR3=0.  5*K*f*H*( )  ♦ 


144 
STR3  -  tertiary  stress 


P  -  density  of  sea  water  (64  lb/ft J) 

H  -  head  of  sea  water  (ft) 

a  -  long  dimension  of  plate  (in) 

S  -  short  dimension  of  plate  (in)  =  distance  between 

longitudinal  stiffeners 
t  -  plate  thickness  (in) 

The  stresses  are  more  severe  at  points  A  and  B  but  we  are 
particularly  interested  in  the  stresses  at  point  B  because  there 
they  will  combine  with  the  stresses  resulting  from  bending  also 
called  primary  stresses. 

Note  that  for  a  panel  under  lateral  load,  the  direction  of 
the  stresses  is,  at  the  edges,  perpendicular  to  them,  so  at  B 
they  are  directed  "longitudinally". 

For  a/S>2  we  have  K=0.685 


Replacing  all  the  numbers  in  the  above  formula,  we  obtain 

S   2 
STR3=0. 152222*H*( ) 

t 
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The  following;  stress  definitions  and  figures  3  and  k   will 
help  understanding;  the  reasons  behind  the  calculations  made  in 
many  of  the  following  sections. 

STRESS  DEFINITIONS 

STH1  -  SHIP  BENDING  STRESS  or  PRIMARY  STRESS 
Since  this  stress  is  directly  proportional  to  distance  from 
the  ship's  neutral  axis,  it  is  substantially  constant  at  any  point 
on  the  cross  section  of  the  bottom  or  top  plating  when  the  vessel 
is  in  an  upright  attitude.  For  distinction,  this  stress  may  be 
termed  a  primary  stress  or  an  area  stress  in  this  location  and 
under  such  a  condition.  While  the  longitudinal  direction  is  imp- 
lied in  the  foregoing  comments,  they  are  equally  applicable 
transversely. 

STR2  -  GIRDER  BENDING  STRESS  or  SECONDARY  STRESS 
Girder  bending  stress  arises  from  the  reaction  of  a  plating 
stiffener  combination  to  a  loading  normal  to  the  plating  such  as 
may  be  due  to  water  pressure.  Following  the  usual  concept,  the 
plating  acts  as  one  flange  of  the  composite  girder  and  so  is 
stressed.  Due  to  shear  lag,  this  secondary  stress  in  the  plating 
is  a  maximum  in  way  of  the  stiffener  and  diminishes  with  distance 
from  it. 

For  the  previous  case  of  bottom  plating  and  vessel  upright, 
the  transverse  variation  of  fore  and  aft  stress  will  be  sensibly 
constant  for  some  distance  from  the  stiffener  but  will  vary 
significantly  through  the  plating  thickness  as  the  plate  thick- 
ness is  a  substantial  part  of  the  distance  to  the  neutral  axis 
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of  the  plate  stiffener  combination.  For  convenience  such  a  stress 
may  be  termed  a  line  stress  although  for  a  wide  spacing  of  stiffe- 
ners  with  consequently  greater  deviation  from  uniform  distribution 
some  inaccuracy  of  concept  is  introduced  by  using  this  term. 

For  girders  continuous  over  several  supports,  such  as  tran- 
sverse bulkheads,  maximum  bending  moments  occur  at  the  supports 
and  will  readily  be  in  the  order  of  twice  the  field  moment. 

STR3  -  PLATE  BENDING  STRESS  -  TERTIARY  STRESS 
Proceeding  in  turn  to  ever  more  elemental  portions  of  the 
ship's  structure  concludes  with  the  simple  plate  panel  supported 
on  its  four  edges.  In  such  a  plate  under  lateral  or  normal  loads, 
bending  stresses  are  created  which  are  variable  from  point  to 
point  on  the  panel  and  through  its  thickness  with  maxima  occurring 
at  one  or  the  other  surface.  Since  any  such  stress  is  effective 
at  a  point  only,  it  may  be  called  a  point  stress  or  tertiary 
stress. 

Each  elemental  panel  of  ship's  plating  is  generally  one  of 
a  repeating  pattern.  Quite  appart  from  any  rotational  restraint 
provided  by  the  panel  supports  consisting  of  floors,  frames, 
longitudinals,  etc.,  there  is  the  restraint  afforded  by  the 
continuity  of  the  plating  and  its  loading.  For  pratical  purposes 
then,  a  condition  of  edge  clamping  exists  with  maximum  stresses 
being  found  at  the  midlength  of  the  longer  panel  edges  and 
secondary  maxima  at  the  midlength  of  the  shorter  edges.  In  both 
cases  these  stresses  are  normal  to  their  respective  edges. 


BENDING     STRESSES    IN    SHIPS 


20 


> 

N.   ) 

~ I  A 

} 

1 

_^/| 

AP 

■ 

1 

FP 

STRESS       DISTRIBUTION       IN 

BOTTOM       SHELL     AT 
TRANSVERSE     WEB         FRAME 


HULL       BENDING      STRESS 


1 


TYPICAL        BOTTOM        STRUCTURE 

OF      A 
LONGITUDINALLY       FRAMED       SHIP 


TRANSVERSE    WEB    FRAME 
OR   BULKHEAD 


a|  °i    °i  °i 


SECTION        A-A 


STIFFENER         BENDING        STRESS 


N* 


i  I  ii  t  iiii  m  itti 


» 

J 
t 

\ 

SECTION        B-B 


PLATE      BENDING      STRESS 

t 


f    I    t"T   I    I    I    I    I    I    ITT 


I  _N» 


FIGURE  3 


21 


From  the  previous  definitions  and  figure  3  we  may  summarize 
in  figure  4  the  effect  of  STR1,  STR2  and  STR3. 
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Fig.  4 
At  point  A 

STR2  is  larger  than  at  point  B 
STR2  is  acting  in  the  longitudinal  direction 
STR3  is  larger  than  at  point  B 
STR3  is  acting  in  the  transverse  direction 
STR1  and  3TR2  are  here  combined  acting  both  in  the 
longitudinal  direction 

At  point  B 

STR2  may  be  acting  in  the  tansverse  direction 
STR3  is  acting  in  the  longitudinal  direction 
STR1  and  STR3  are  here  combined  acting  both  in  the 
longitudinal  direction. 
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Usually  STR2  is  attributed  a  constant  value  and  can  be 
taken  into  account  by  automatically  decreasing  the  design  stress 
by  the  value  assigned  to  STR2. 
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THE  BENDING  MOMENT  EQUATION 
The  bending  moment  equation  used  is  derived  from 
BM  =  A*AL/35      ton*ft 


knowing  that 


and  that 


and  using 


we  obtain 


A  =  V/35    ton 
V  =  B*T*AL*Cb 
Ch  =  0.75 


0.75 

BM  =  *B*T*AL*AL     ton* ft 

35*35 


where  AL  =  is  the  length  in  ft 

B  =  is  the  beam  in  ft 

T  =  is  the  draft  in  ft 

C^  =  is  the  block  coefficient 

or  multiplying  by  2240  we  obtain 

BM  =  1.37143*B*T*AL*AL      lb*ft 

Note: 

V  is  the  displacement  volume  and  £±    is  the  corresponding 
displacement  weight. 
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MOMENT  ACTING  IN  ONE  DIRECTION 

NO  RESTRICTION  IMPOSED  ON  THE  COMPRESSION  STRESS 

Considering  the  box  girder  subjected  to  a  bending  moment 
acting  in  hogging  BM,  we  will  observe  as  a  consequence  that  the 
top  will  be  subjected  to  tension  and  the  bottom  to  compression. 
These  stresses  are  given  by  the  already  introduced  formula 


BM*Y 
STR= 


or  in  this  case,  being  interested  on  the  top  and  bottom, 
BM*DIBTOP 


STH= 


SBH= 


I 

BM*DIBOT 


In  this  introductory  study  it  was  decided  not  to  limit  the 
value  of  the  stress  in  compression.  Only  the  value  of  the  stress 
in  tension,  here  STH  will  be  requested  to  be  equal  to  the  value 
chosen  for  the  design  stress  in  tension  (or  smaller). 

SBH  will  be  a  consequence  of  the  value  of  DST  and  also 
of  the  S  value  chosen. 

So,  we  will  basically  have  3  equations 

DST  =  STH«*TOP    Bending  equation 

S  ■**  SBH  <*♦•  BOT     Bending  equation  +  Bryan  formula 


STH 


BM 


( 


SBH 


COMPRESSION 


TENSION 


Fig.  5 
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The  three  equations  that  can  be  obtained  will  allow  one  to 
determine  the  three  unknowns  which  are  BOT,SBH,TOP. 

By  replacing  in  the  equations  already  mentioned  the  expre- 
ssions derived  previously  for  DIBOT,  DITOP  and  I  and  solving  for 
BOT,TOP  and  SBH  we  obtain  after  a  lengthy  but  simple  process  the 
equations  in  their  final  format  as  we  will  see. 

Note  here  that  in  order  to  solve  the  three  equations, 
the  thickness  of  the  side  plates  TS  was  made  equal  to  the  larger 
TOP  or  BOT;  since  TOP  and  BOT  were  not  known  and  had  to  be  deter- 
mined after  deciding  on  TS,  the  equations  had  to  be  solved  twice. 
First  the  equations  are  solved  assuming  that  TOP>  BOT  this 
simply  implying  that  the  neutral  axis  is  going  to  be  chosen  alrea- 
dy closer  to  the  top  and  TS  =  TOP.  Second,  the  equations  are 
solved  assuming  BOT > TOP  thus  making  TS  =  BOT. 

As  will  be  seen,  only  one  solution  is  possible,  the  other 
being  rejected  by  giving  a  result  that  does  not  agreee  with  the 
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assumption  just  mentioned. 

TOP  ,  EOT  ,  S   (inches) 

B    ,    D    ,    DITOP  ,  DIBOT  ,  AL   (ft) 

DST    ,    STH    ,    SBH    ,    E      (psi) 

BM     (ton*ft) 

I      (in*ft)2 


560*BM* ( B*BOT+D*TS ) 

DST£  STH   = 5- 

D* ( 3*B*B*TOP*BOT+2*B*D*TS* ( TOP+BOT ) = ( D*TS )     ) 

560*BM*(B*TOP+D*TS) 

SBH   =  


Same  denominator 

BOT  =  \J  G  >J  SBH 

(S)2*(l-(POI)2) 


G  =  0.30396* 


DST    B*BOT+D*TS 


E 


SBH        B*TOP+D*TS 

Making 

9*B*B*D 
A   = 


C   = 


140* 

►BM 

3*B*D*D 

70* 

►BM 

3*D*D*D 

F  =• 

140*BM 


we   will   obtain,    by  assuming  first  T0P>  BOT 
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D*(B+D)*TOP 


DIBOT  =• 


(2*D+B)*T0P+B*B0T 
D*(B*BOT+D*TOP) 


DITOP 


Same  denominator 

(  2*B*D+D*D )  *TOP*TOP+  (  2*B*D+3*B*B )  *TOP*BOT 


I  =  4*D*D*- 


Same  denominator 
B* ( BOT ) 3 


TOP 


BOT  = 


( B+D ) *G*DST-D* ( BOT ) 2 
TOP* ( 1 2*D- ( C+F ) *DST*TOP ) 
DST*(A+C)*TOP-12*B 


Before  going  to  the  equations  obtained  when  BOT >  TOP 
is  assumed,  we  call  the  attention  of  the  reader  to  the  fact  that 
the  values  of  BOT  and  TOP  have  to  be  obtained  from  the  two  above 

equations.  Substitution  of  one  into  the  other  would  result  in 

fch 

an  equation  of  degree  higher  than  at  least  5   degree,  unsuitable 

for  direct  solution.  As  will  be  later  exemplified,  the  solution 
is  here  obtained  by  assigning  values  to  either  one  of  BOT  and  TOP 
and  through  the  two  equations,  we  obtain  the  respective  values 
for  TOP  and  BOT.  The  ploting  of  these  values  BOT  vs  TOP  will 
give  one  line  for  each  of  the  two  equations.  The  intersection  of 
the  tv/o  lines  thus  obtained  allows  the  solution  to  be  read  from 
the  plot. 

This  solution  will  be  acceptable  if  the  intersection  lies 
in  the  sector  TOP>  BOT  otherwise  the  intersection  must  be 
rejected. 

Assuming  now  BOT >  TOP  we  make  TS  =  BOT  and  obtain: 


2b 


D*(B»TOP+D*BOT) 
DIBOT  = 


(2*D+B)*BOT+B*TOP 
D*(B+D)*BOT 


DITOP  = 


Same  denominator 

( 2*B*D+D*D ) * ( BOT ) 2+ ( 2*B*D+3*B*B ) *TOP*BOT 


I  =  l±*D*D*- 


TOP  =• 


BOT  =■ 


Same  denominator 
( B+D ) * ( BOT ) 3»D*G*DST*BOT 

B*G*DST 
12* ( B+D )  - ( A+C ) *DST*TOP 

( C+F )  *DST 


In  this  form  the  equations  can  be  used  in  a  very  simple 
computer  program  that  will  give  the  set  of  values  TOP  vs  BOT 
suitable  for  ploting.  At  this  stage  we  are  hoping  that  if  the 
intersection  obtained  in  the  situation  where  we  assumed  TOP>  BOT 
is  acceptable , the  two  lines  resulting  from  the  case  BOT>  TOP  will 
not  intersect  or  if  they  do  that  the  intersection  must  be  rejected 
due  to  the  fact  that  they  will  violate  the  assumption  from  which 
they   were  derived 

As  can  be  observed  from  the  several  cases  run  in  the  compu- 
ter, the  case  where  two  acceptable  solutions  exist  never  happe- 
ned and  it  can  only  happen  if  TOP  =  BOT  where  then  both  solutions 
are  acceptable  and  coincide  as  in  the  case  of  one  of  the  follo- 
wing runs. 

In  the  following  pages  the  computer  printouts  are  included. 
The  first  run  is  explained  in  detail  to  make  the  interested  rea- 
der familiar  with  the  problems  encountered  in  determining  the 
intersections  and  the  respective  thicknesses. 
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The  other  runs  are  here  presented  as  well  to  serve  as 
proof  that  a  simple  solution  can  always  be  obtained  and  also  to 
give  an  idea  of  the  range  of  thicknesses  expected  in  each  case. 

In  all  runs  only  a  bending  moment  acting  in  hogging  was 
considered. 

The  program  used  is  a  direct  application  of  the  equations 
here  derived  to  a  particular  set  of  values  chosen  as  the  input. 
From  the  results  obtained,  some  conclusions  were  drawn  that  a- 
llowed  for  the  improvements  made  in  following  sections. 

No  further  explanation  on  the  program  shall  here  be  given 
because  it  is  self  explanatory  in  the  comment  statements 
included. 
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Prom  these  lists  of  thickness  values  printed  by  the  com- 
puter, we  can  see  that  the  acceptable  ones  are  those  marked 
with  a    sign.   Plotting  these  corresponding  pairs  of  values 
properly,  as  is  next  exemplified,  allows  for  one  to  read  from 
the  intersection  of  the  two  lines,  which  falls  within  the  range 
set  by  the  assumption  that  originated  it,  the  values  of  thick- 
nesses at  the  top  and  at  the  bottom  which  satisfy  the  model 
we  started  with. 
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o.%  ■ 


on-- 


0,6.. 


top 

thickness  (inches) 


0.2>  ■• 


0.2. 


0-1 


-4 1- 


-I u 


0.1         o.l 


0.3      o-4      °=> 

bottom  thickness ( inches ) 


We  can  see  that  the  thickness  values  are  easily  read  from 
the  intersection  obtained  and  they  are  acceptable  because  they 
fall  within  the  range  assumed  to  derive  the  formulas  that  ori- 
ginated them  namely  TOP>BOT. 

Next  will  be  plotted  the  case  BOT>TOP  and  as  can  be  seen, 
again  an  intersection  is  obtained  which  in  this  case  is  not 
acceptable  because  it  occurs  within  the  range  TOP>  BOT  contrary 
to  the  assumed  BOT > TOP  that  originated  it. 

With  a  little  more  experience,  this  can  be  seen  immediatly 
from  the  computer  printouts,  thus  avoiding  loss  of  time  doing 
a  plot  that  will  lead  to  a  nonacceptable  result. 
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top  thickness    (inches) 


0-1  .. 


0-X 


0,3  0.«f  C\S~  (XG 

bottom   thickness    (in.) 
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Now  having  obtained  the  values  for  the  thicknesses  from 
the  plots,  we  verify  them  through  the  use  of  the  next  program 
based  on  the  formulation  previously  derived.  As  can  be  noted, 
we  obtain  an  acceptable  value  for  the  stress  at  the  top  (since 
we  are  considering  a  hogging  bending  moment)  which  we  said 
before  in  the  derivation  of  the  equations  should  not  exceed  the 
design  stress  in  tension.  In  this  case  the  design  stress  was 
DST=20000  psi  and  we  see  that  the  stress  at  the  top  STH=20037.5 
the  37.5  coming  possibly  from  inaccuracy  when  reading  the  thi- 
cknesses from  the  plot,  and  computer  roundoffs. 

Note  that  the  stress  at  the  bottom,  since  there  was  no 
limit  placed  on  it  except  the  Bryan  formula  to  resist  buckling 
caused  by  compressive  inplane  loading,  is  attaining  a  too  large 
value  which  is  unacceptable  and  unrealistic.  Here  again  there 
should  be  no  difference  between  SBH  and  SBCRIT  and  the  difference 
that  can  be  noted  is  also  due  to  inaccuracy  in  reading  the  thi- 
cknesses from  the  plot  and  computer  round-off s. 
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SUMMARY  FOH  THE  PRECEDING  CASE 


Length 

= 

500  ft 

Beam 

= 

AL/8 

Draft 

= 

B/3 

Depth 

= 

AL/9 

DST 

= 

20000  psi 

S 

= 

10  inches 

Material 

- 

STEEL 

TOP 

= 

0.468  inches 

BOT 

= 

0. 166  inches 

STH 

= 

20037.45 

SBH 

= 

30434.89 

SBCRIT 

= 

29886.80 

DITOP 

= 

22.06  ft 

DIBOT 

= 

33.50  ft 

We  can  already  notice  here  that  the  stress  in  the  bottom 
is  attaining  a  too  large  value  due  to  the  fact  that,  since  no 
limit  was  imposed  on  it,  it  is  dependent  only  on  the  chosen 
value  for  S  and  the  respective  capacity  to  resist  buckling  pro- 
vided by  S  and  BOT  while  at  the  same  time  stressing  the  top  plate 
to  its  maximum  allowed  DST. 

In  the  following  two  cases  only  the  spacing  S  between  the 

longitudinal  stiffening  members  will  be  allowed  to  change  and 

the  results  obtained  are  simply  stated  as  previously  exemplified 

and  were  determined  in  a  similar  way. 

From  the  results  obtained  we  can  note  the  effect  caused 
by  the  change  of  S. 
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**  < 

2  2  6  3 

~z-?  = 

C  .  c  7  "  7 

3  c~  2 = 

0  - 

"923 

'CD2  = 

? .  5  5  7  z 

0,  n  -  p  _ 

7  ■ 

?  •.  a  7 

—  r  = 

P  .  6  0  7  0 

3C~2  = 

-7 « 

TOIf 

TnppS 

5»  •  f  55"f 

30"?  = 

-pi 

1  "4  h 

TPPps 

r>  ,  7  r/ 1?  1/ 

RO-?  = 

•pi. 

4  i»87 

TrPp  = 

V  .  7  ~  7  '/ 

RO-Pa 

-1?. 

1«»8 

TOD?s 

7  .  E7PI7 

Rn*?  = 

-71 « 

p?AI? 

T  r  ■  p  p  _ 

0,8500 

30'2  = 

-o  . 

2^12 

T  o  P  g  « 

1  .97i^P! 

BP~2  = 

-r/. 

t>  ^50 

T  ■;  D  2  = 

7  .  °570 

30~2- 

-01 

3276 

TQpgB 

1  .0000 

B0'2- 

»0. 

3693 

^5 


top   thickness    (in.) 

O.fe 

\ 

1         / 

/ 

0.>'. 

\     / 

/   / 

0.4. 

/'  j 

0.3 

/'         I 
/               j 
/                   j 

0.1 

/                       f 
/                            I 

/'                          J 

0-1  ■ 

1 1 < V 

1 »- 

1 

O.i       02.      0-i      0.+       0.b'      0.(0 
bottom  thickness  (in.) 
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:nd 


h-tto- 

ASSi  MFD 

LARbFR  THAN  TOp 

T  )P3« 

•  1 

.777?. 

B  0  T  3  ■ 

n 

.  11-10  0 

T->P  ,= 

-e 

.;  94  9 

BOT  3  ■ 

0 

p  }  f  <  7  0 

T^P3. 

-. 

.  0  K  f,  7 

B0T3- 

0 

.  p  . '  0  0 

T  -;  p  3  ■ 

-■/ 

.'/'Ht 

BCT3«= 

»?500 

T^Pje 

z 

.  .'V, "6 

BO-3* 

/■» 

•  3000 

T  )P3« 

■ 

•  1769 

50~3  = 

r> 

►  3500 

TCP3- 

r 

•  37P9 

^Q"3^ 

7 

•  4  2  ?  0 

T  Q  u  7  3 

-■ 

.<-37i 

B0'3» 

0 

'450  0 

T  c  P  3  ■ 

/ 

.  q  7  g  2 

BC~3  = 

0 

•  50  7  0 

T"P-;e 

l 

•  4047 

5CT3  = 

7 

•  5  5  0  0 

T  ^  p  3  ■ 

l 

.^25^ 

RT*3e 

0 

-6-77 

T (  c-je 

g 

,5478 

=^T3  = 

7 

-6CZ7 

-  -:  ■ 

3 

»?816 

^CT3  = 

7 

'77C? 

T  '•  t  "5  — 

*» 

1349 

30T3" 

7 

>7^00 

T  '  '-'"5  x 

K 

1163 

B0T3« 

0 

8000 

T  H  D  T  = 

A  , 

?3fc0 

BO"?  = 

2  * 

8  c>  0  ? 

T   :"  Os 

7  . 

4970 

30T3" 

7< 

9?  £7: 

''Tis 

w  , 

^137 

^CT3  = 

7  ' 

0^00 

'  '  "  3" 

10, 

49?5 

B0T3«« 

1  ' 

0  |  0  7 

T  -  If- 

(7.  , 

1000 

B0T4« 

?< 

on  fcl 

T  j  P  4  ■ 

7  i 

1  5(  .' 

5?"^  = 

0  ' 

97  31 

''Ol,: 

7  . 

P000 

B0T4« 

7' 

P101 

T  '  -  4  = 

'/. 

2500 

3£T4a 

?' 

71711 

T  P4s 

?  . 

3  2  0  0 

R^T  4  s 

7< 

6240 

TnP;  = 

e. 

3500 

POT4_ 

e  • 

5310 

*   Pfc« 

■V  i 

4.000 

pn't  = 

0i 

4379 

T'.Wi»  = 

. 

4  5/? 

B0"4  = 

0i 

34  4  9 

'  P  4  ■ 

5000 

g  q  t  u  m 

7' 

?519 

T^u 

7' . 

R500 

RnT/.= 

0« 

15«9 

T^pii  = 

6000 

qo^4_ 

7' 

P!f58 

T  r  p  /» = 

e . 

6  5  7  0 

yn-4  = 

-0, 

7?7H 

TrP4  = 

f . 

7000 

RO-4K 

-7  < 

1P02 

T  0  P  4  a 

~ , 

75  7'/ 

boT4  = 

.0, 

PI  3  3 

Tr?  D4  = 

7  . 

R  0  0  3 

ROTte 

-0. 

3?63 

T  r  p  4  = 

?  . 

8500 

50-4  = 

-7" 

?  a  P3 

ToP^a 

7  . 

900  0 

B0T** 

-7  < 

4r  ?3 

TrOfc- 

0i 

•?50>2 

9n-4  = 

-?• 

5851 

Tc^" 

1  • 

0000 

B0T4« 

-^ . 

6784 

OPERATING  SYSTEM  -ErSION  1  REVISION   ill?  03/04/7*  GENERATED   0-5/06/74 


JOB 

173 
67 
00 
20 


-ANGLING  CHARGE    s   ,35 
LINES  PRI\TEn  RR?  *  1.25 


CARDS  REA- 
c  L  0  T  T  E  R  V '.  C T T R S 
- CGEl  70  SECONDS 
^OOEL  BZ  SECON!  3 


*  1  .c 

*  .85 
•?5.00 

•12.50 

TOTAL 


K  LN 

<  CO 

1  ;-  7  >' 

h  -  u  n 


CMA^ 


•35 

•  22 

•  IP 

•  0? 

•  14 
.07 

•  81 


(El-  4oj  lW3i   LOGGED  TUT  75/06/7*  15120.   *   13.  ?9  LtET  AFTER 


18  LOGINS. 


top  thickness    (in.  ) 
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or 


O.b- 


0.b 


o.«*. 


0.3. 


o-Z.. 


0.1 


-I u 


0.1 


0.2  0.3  o.M-  o.d  o.6 

bottom   thickness    (in.) 


strenght   Or   ships-on   the   nfut^al   axis  location 
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C 

- 

18 

c 

24 

•.' 

c 

- 

3: 

35 

c 

3? 

- 

s 

r 

- 

- 

74 

- 
- 

■ 

I  - 


•  U 

:  •  : 

DEN 


-AL*  IS  THfc  LENGHT 

al  c5aa . 

-••  =  A  L  /  8  . 

,  DRA*  IS  ThE  draft 
DRA"B/3t 

[)=AL/9« 

*PQI»  IS  THE  POISSON'S  RATIO  OP  THE  MAtfRIAl 

FOR  STEEL 

PQI"Q«3 

-E.»  IS  THE  YOUNG'S  MODULUS  OF  THE  MATERIAL 

steel 

"  a  353  ■»ltti**6 

..:-t*  is  The  design  stress 


s  • 

* 


IS  THE  SPACING  OF  LONGITUDINAL  STIFFENING  MEMBERS 


r'    _     —    -n 

i»3« 


,  '-    ••   : 


IS 

«s ( e  * i r 

i»?  •  3C  3  9 

«  7  Q  P  »  '; 

-     '  T -     IS 

T    P«s0.4] 

Q0T=?»41 

•  =  w,  -> .  *  ;- 

'•=--( 3 

«  c  T H *     IS 

STh*A/DE 

-.  i  TE  •  5  # 

6{)    ■  OPM AT  (  • 

■  a  -  <-  ?  .  »  :- 

-  -''-•*     t  S 

•'•  J=C/DE 

'  TT=-f5» 

96     r-r'AT  (   I 
•S3CRIT* 

-•..pOJT.-. 

r I T{  '  5  / 

116    FOR M ft T (  • 
peJ.*n*T 
- jtoP. 

•-,  J  T  P  D  B  '     ,< 

"  T  -   f  r:  , 
t  (      FOP H  A  T  (  I 

DI80T«D* 

- IT? ( 5j 
50    POR«AT ( i 

I  ND 

c  3?* [v] 

C  3*CCV3 

I  37C  Cv] 
fi3H^  CV] 


TH£  3 E N D I  '-■■  G  m  Q  M  r  N  T 

A  -  \  L  *  A  _  •"/  t  7  5  )  /  (  3  5  .  *  3  5  .  ) 

6«<3*S*(1»-P0I*P0I)/E 

T-E  THlCKNirSS  OF"  THE  TOP  PLATING 

the  thickness  of  the  bottom  Plating 
25 

25 

m*(9#83T+D*T0P) 

.«-.  «-«TO"*"DT  +  ^.*3>«O»'TOp«>(T0P  +  HOT)  +  (D*D*TO»*T0P)) 

T  -E  STRESS  0!"  TOP  DUE  TO  HOGGING 
N 
60 |STH 

'  t  20X ,  '  S  T  LJ  n  '  F  1  7  .  3  ) 
fl  »  |  B  *  T  3  P  + [J  *  T  3  •-'  ) 

*..r  5^rq;  qn  gQT  DUE  TO  HOGGING 
N 
90 |SBH 

1  >20X,  •S3'-'«'  P10»3) 

IS  T-JE  CRITICAL  STRESS  ON  BOT 

n  T  ,  r,  ,,  t  /  3 
1  -   )  S  ?  C  R  I  T 

i  ,  i  7  x  ,  '  s  --  3  R  !  T  =  '  r  1  C  •  3  ) 
0P+8*(T0P+acT) 
IS  THE  DISTANCE  OF  Tl-IF  TQP  FROM  THE  NEUTRAL  AXIS 

( B*B0T*D*T0P > /P 

u  ;■.  )''TT"n 

'  ,  J  R  X  ,  '3XTOD*'F10*iH 

IS  TwF  DISTANCE  OF  THE  BOTTOM  FPOM  THE  NEUTRAL  AXIS 

(B*ToP+D*TO° )/P 
50 ) 0 I BOT 
•jlRXj'DIBOT-'FlfltiH 


AL 

t 

G 


0  3  2  C  C  V  ]  B 

'/;•'-■    [S3  .  R 

7^>rv;  tqp 

?  1  7  <♦  C  L  ]  60 


K334[V]  DRA 

035CIV]  CST 

0  3  9  0  C  V  J  e  0  T 

0  0  0  0  C  S ]  9  I 


033CCV]  D 

036^CV]  S 

0  3  9  4  [  V  ]  A 

03B3CV3  C 


'•TkE^CHT     Of"     SHIPS-PN     ThE     NMjTi-AL     AXIS     LOCATION 


ctv?    SRM 

cv:   cito"3 
>Fc; 


:-lE?CL] 
P2AS CL3 


90 


■/3cv  cv] 
03CCCV] 


SRCR'T 
DIt.<0T 


0228 CL]     110 
0304 [L]     50 


03C4 CV] 
00'/ir/  CS] 


^9 
PAGE 

i 

V 


cgr     els: 


•MAIN* 

2  6  *  0 

.V 

2FQr 

•  COMP 

2A6E 

(u  I 

2440 

•  R 

29b8 

.Zi 

«6 

2480 

.  A 

29"4 

»MES 

2832 

.  W 

263A 

EXP 

27C0 

Air 

?b32  ALOG 

28F8 

•  PAR  - 

2954 

•  5 

?63A 

AE.XP 

2950 

•ERCNT 

29?A 

*8 

?95C  »0 

2A32 

•  U 

3P-  i 

-points  : 

.R 

?4- 

.  A 

253? 

A  LOG 

?a?a 

ALX^ 

263A 

EXP 

2700 

Air 

••■ 

2*90 

•  COMP 

28B6 

$6 

2  ^ES 

i  R  A  R  G 

292A 

*8 

2954 

•  5 

•  ' 

295C 

.  ERC  T 

295F 

.0 

29B4 

•  NFS 

2A3E 

.u 

2  AC,  4 

•  V 

H 

• ■  locks: 

E 

D  SUBROUTINES! 


,>.qr:u  ADORLSS  2?70 
BEGINS; 


STH=  2  1037. 113 

SBH«  2T37.113 

SpCRIT=  2  5e5.348 

DlTOP=  27*7778 

OlBCT=  27*7778 


IF  VIC  OPEkATI-  G  SYSTEM  VERSION  1  REVISION   012  03/04/74  GENERATED   05/06/7<f 


HANDLING  ChAPGE  *   .3b  /  JOB  .35 

LINES  PPI-  TED  PR2  *  1.25  /  K  L.N  «11 

5*  CA9DS  READ  £  1,56  /  <  CO  »0.8 

?0  PLCTTtrP  VrCTORS  *   .25  /  1'  30  .00 

13  *ODEL  7/    SECONDS  »25.0Z  /  HOUR  .09 

"EL  80  SECONDS  M  2  .  5  •'  /  HOUR  .00 

TOTAL  CHARGE  *     .63 


EIR  k3?     I473t   LOGGED  "Ut  05/06/74  15:46.   *   11. lft  LEFT  AFTER 


21  LOGINS. 
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SUMMARY 

Length 

= 

500  ft 

Beam 

= 

AL/8 

Draft 

= 

B/3 

Depth 

= 

AL/9 

DST 

= 

20000  psi 

S 

= 

30" 

Material 

- 

STEEL 

TOP 

= 

0.4125" 

BOT 

= 

0.4125" 

STH 

= 

20037.11 

SBH 

= 

20037.11 

SBCRIT 

= 

20505.43 

DITOP 

= 

27.78' 

DIBOT 

= 

27.78' 

Note  that  this  is  the  case  where  two  acceptable  solutions 
exist  although  they  end  up  being  the  same  because  TOP=BOT. 

Note  that  the  neutral  axis  is  located  half  way  between 
top  and  bottom  and  the  stress  is  the  same  at  the  top  and  at  the 
bottom,  as  it  was  expected  beforehand. 


STRENGHT  Of  ships-on  the  neutral  axis  location 
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c 

c 

r 

c 


*A|_*  IS  THE  lENGhT 

fcL«500« 

B  =  AL /8 • 

*DRA*  IS  ThE  DRAFT 

D  R  A  ■  B  /  3  • 

D  =  AL/C|. 

*P0I*  IS  ThE  POISSO.M'S  ratio  of  the  material 

FOR  STEEL 
pOI*0«3 

*E»  IS  TmE  YOUNG'S  MODULUS  OF  THE  MATERIAL 

FDR  STEEL 

£  =  3 '/  .  »  1  /  .  *  *  6 

•DST»  IS  THE  DESIGN  STRESS 

0  S  T  =  S  v- :' .  /  • 

*5«  IS  THE  SPACING  OF  LONGITUDINAL  STIFFENING  MEMBERS 

*SM»  IS  th:  bending  moment 

5M=(B»,DRA*AL*AL*««75)/<  35 .  «  35  ■  I 
2=0 • 30 396 "S * S * ( 1 • -PO I  * PO I ) /E 

A  =  a  .  *  r-  ■  B  *  D  /  (  l  k  7  •  *  BM  ) 

0  -  3  '  •  r*  -  :-i  *  0  /  (  7  Z  •  *  B  v  ) 

F=3«*D*D*D/(l<f,  •  *  B M  ) 

THICKNESS  Cr  TOP  ASSUMED  LARGER  THAN  BOTTOM  THICKNESS 

10  FORMAT (» 1 » ,25X, 'TOP  ASSUMED  LARGER  THAN  BOTTOM') 
.JrP*     TS  TwE  THICKNESS  OF  THP  TOP  PLaTj.\G 

•°CT*  IS  THE  THICKNESS  OF  THE  BOTTOM  PlaTING 
BOT!  =»1  • 

11  ~  <~  o  <  - b * a 0 T 1  *B0Ti*30Ti/<  (B+D)*G*DST»D*BoTl*BOTi  ) 
W  F :  I  t  r  (  5  ,  p ',""  |Topl^POTl 

20  FCRMATM  %  25X,  «  Tr>Pl  =  »  F9  »4*  5X#  '  SOT1*  ♦  Fg,  <t  ) 
aOT!  =nOTi  *:■  #  ?5 
•  r  ,  n n t  1  .  r, T  .  2  .  )  G 0  TO  3 0 
G  C  "0  II 

30  TCP?=0^02 

31  nOT?=TOPP*r 12»«D-(C+F) *DST*T0P2)/(DST*( A+C) *T0P2-12«*B) 
W  R I T  r  I  5  /  * '«""  '  TO  p  P  /  n  0  T  ? 

%e  FORMA"T(«  ',  ?^x,   '  T0P2«  .  pg,  Ly,  5X,  •  BOTP  =  ■  Fg.  4  ) 

TnP2oTOP?+O.02 

jP( T0P2-GT . 1 . >G0  To  5T 

OC  "r  31 
c.v    C5KTIMJE 

THICKNESS  ^F  R^TTnv  ASSUMED  LARGER  THAN  TOP 

HRlTE<5/60) 

60  pORfATC  •  1  •  rP.S^t  'BOTTOM  ASSUMED  LARGER  THAN  TOP') 
=  CT">  =  1  • 

61  T 0P3« (  lt:+C)*nPT3*H0T3*30T3-D*G*0ST*B0T3)/(B*G*DST) 
&.RIT  '  B,  7?  )ToP3/BOT3 

70  rORMAT('   •  ,  2--x,  •  Trp?=i  F  9  •  <♦  >  5  X  ,  •  B0T3=  ■  Fg,4  ) 

.-.  c  t  3  ■  f.  p  t  9  ♦  - .  7  5 

TF( B0T3.GT.2, (GO  TO  80 

GO  TO  61 
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•402 

H<fE 
>7? 
'49A 

5*BC 

.n[v] 

J?tV] 

•   cu 

SdCU 
-    CV3 

7c  cv: 
xeo 

DP 


80 
81 

90 


T0P4=fc .02 

nOT4»{  12.*{B+D)-(  A  +  C)1[DST*T0P4 

WRITE (5^98 |T0P4jB0T4 

FORMAT( •      * t2*Xt  • TCPts 'F9.<*/5X, 

TOP4»TOP4+0i  Z>S 

IF(T0P4»GT»1.  )GO    TO    1(10 


)/ ( ( OF ) *DST ) 

'BOT<*=  ■  F  S.h  ) 


GO    TO    81 
100    CONTINUE 

END 
.  U  P«*»C 0  CVi     & L 

POI  04E8CVJ     E 

B  M  0  5 1 C  [  v  ]     G 

\Z  0000 CS]     «I 

20  0  ?  2  c  c  l  ]    30 

A0  ^?Fc  r.L]     50 

T  0  P  3  0  3  p  ?  r.  L  3     70 

P0T4  7A72LL2     9^ 


506RAN   labels: 

2071     *MAir  * 

Z6DE    ALOG 
cP0^     .0 


ntry-pcints  : 

?5EC  »R 

?9Dr  «W 

2B0*  .ZERO 

?C1A  f»I 


OfMC^-RLOC<S: 


BC0C  «V 

262C  »A 

2  A  A  <<  .  R  A  R  G 

2RDE  »U 


?6?C     »A 
2a3C    »COMP 

2  P  0  8     .ERC\T 


DEFINED    SL'BROljTIN'ES: 


SFER    ADnPESS    2^7? 

ecutio.  begins: 


04C8EV] 

0  0 00  CS ] 

0528  rvl 
053C  [V] 
055<i  CV3 
0312 CL- 
0402 CL3 
f+PC  CL- 


B 
•  R 

A 

B0T1 

T0P2 

6'/ 

80 

1  0  2 


04D0  CV] 
04F8  CV] 
0530 CV] 
0 1 7  6  C  U  3 
0234  EL] 
056CCV; 
0  5  7  4  [  V  l 
0  0  0  0  [  S  ] 


DRA 
DST 
C 
11 

31 

B0T3 
TOPi» 
.V 


2A3C  «COMP 
2B60  iMES 
2  3  0  g  •  5 

3036 


26DE  ALOG 
2  A  6  2  *  ft 

2B0A  .0 


2C1A  PI 
29DE  «^ 
27E6  AEXP 


27E6  AEXP 
2AA4  . R  A  R  G 

2660  .MFS 


04D8  CV] 
0500 CV] 
0  5  3  8  C  V  ] 
0  5  4  0  C  V  ] 

055C  CV] 

0  3  k « :  l  ] 

0*t0A  CL3 


25EC  .R 
27E6  EXP 
2B08  .ERCNT 


27E6  EXP 
2AD6  $8 
2BEA  .U 


D 

S 

F 

TOPI 

BOT? 

61 

81 


2B0  4  .Zl 
296C  AT 

2AC6  $8 


296C 

2B00 
2C10 


AI' 
.5 
.  V 
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TCP  A 
T  -  P  1  = 
T  0  P  i  ■ 
TQPi  = 
T(  Pje 
TcPi  = 

T0Pl  = 
TnPi  = 

T  q  p  i  = 
T  o  P  i  ■ 

T0P1  = 
T^Pi  = 
TC.Pi  = 
T0Pl  = 
TrP<  = 
TcPl  = 
T  r* p  I  = 

ToPi* 

T  r  P  i  s 

TrPi« 
TrP;  = 

T  C  p  1  c 
TrPps 

Trpg  = 
Tr  p^  = 
TrPpr 
T«Dp  = 

T0P2  = 

T"Pj)S 
T^Pp; 

TpP?  = 

T0pp_. 

T  c  P  g  « ■ 

Trp?  = 

TCPr  = 
TC.PP  = 

TrPpa 

TCP2  = 

TCp?  = 

Trpp  = 
TqP2  = 

T0p2  = 
TPP2  = 
TCp2  = 
TCP2  = 
TCP2  = 
T^Pp  = 

Tcp?«= 

T:P?  = 
t  cp?  = 

TnPpB 

TCP?  = 
TCP2  = 
TCP2  = 

Top?« 

TrDpc 


SSUMED  LARGER  THAN 

0.?423  BCT1  = 

0t?868  B0T1» 

0.3378  BCT1= 

0»3960  B0T1- 
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Note  that  in  this  case  by  simple  inspection  of  the 
coordinates  of  the  points  to  be  ploted  we  can  see  that  the 
intersection  will  occur  outside  the  acceptable  range. 

Note  also  that  in  order  to  obtain  the  coordinates 
of  the  points  within  a  range  of  interest  where  the  intersection 
may  occur,  different  ranges  and  increments  are  considered  in 
each  program.  This  could  be  improved  if  at  that  time  these 
printouts  were  available  and  a  better  idea  of  the  thicknesses 
that  should  be  expected  existed. 
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SUMMARY 

Length 

= 

500' 

Beam 

= 

AL/8 

Draft 

= 

B/3 

Depth 

= 

AL/9 

DST 

= 

20000 

psi 

S 

= 

120" 

Material 

- 

STEEL 

TOP 

= 

0.03" 

BOT 

= 

1.1V 

STH 

= 

19771. 

88 

SBH 

= 

9579. 

87 

SBCRIT 

= 

9788. 

37 

DITOP 

= 

37. ^2' 

DIBOT 

= 

18.13' 

Note  in  this  case  the  large  value  chosen  for  S  and 
associated  with  it  the  large  value  required  for  the  thickness 
at  the  bottom  to  avoid  buckling  which  nevertheless  does  not  re- 
sist stresses  larger  than  9788  psi  in  compression. 
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If  we  look  now  at  these  three  cases,  all  for  the  same 
material  and  set  of  dimensions  (same  bending  moment)  we  can 
see  how  the  thicknesses  and  neutral  axis  position  changed  as 
the  value  of  S  was  changed. 


CASE 

s 
inches 

TOP 
inches 

BOT 
inches 

SBH 
psi 

SBCHIT 
psi 

DITOP 
ft 

DIBOT 
ft 

1 

10 

0.468 

0.166 

30000 

30000 

22.06 

33.05 

2 

30 

0.4125 

0.4125 

20000 

20000 

27.78 

27.78 

3 

120 

0.03 

1.14 

9700 

9700 

37.43 

18.13 

In  this  table,  we  can  more  readily  make  conclusions  on 
the  influence  of  the  value  chosen  for  S  on  the  neutral  axis 
location,  and  also  a  quick  look  at  the  resulting  values  for 
the  thicknesses  may  give  an  idea  of  the  weight  associated  with 
each  case. 


In  the  following  cases  the  dimension  parameters  are 
changed  (increased)  and  consequently  the  bending  moment  and  the 
thicknesses  of  plating  will  be  higher. 


BTREnghT  Of  ships-on  the  neutwal  axis  location 
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SUMMARY 

Length 

= 

1000  ft 

Beam 

= 

al/5.75 

Draft 

= 

B/3.3 

Depth 

= 

AL/14 

DST 

= 

18000  psi 

S 

= 

30" 

Material 

- 

STEEL 

TOP 

= 

5.12" 

BOT 

= 

0.62" 

STH 

= 

18085.02 

SBH 

= 

W73.71 

SBCRIT 

= 

46323.82 

DITOP 

= 

19.55' 

DIBOT 

= 

51.87' 

We  might  say  here  that  the  value  chosen  for  S  is  small 
and  this  accounts  for  the  compression  at  the  bottom  to  reach 
such  a  high  value. 


In  the  next  case,  the  value  chosen  for  S  is  increased  to 
60  inches. 
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Since  sometimes  it  is  desirable  to  use  aluminum  alloys 
in  the  structure  of  ships,  due  to  their  light  weight,  high 
corrosion  resistance  and  non-magnetic  characteristics,  in  the 
following  runs,  the  material  considered  is  aluminum.  The 
necessary  change  is  made  in  the  input  to  account  for  the  pro- 
perties of  this  material. 

When  considering  aluminum  as  the  structural  material, 
we  must  keep  in  mind  that  the  aluminum  hull  will  be  more  flexi- 
ble than  a  corresponding  steel  ship.  The  modulus  of  elasticity 
of  aluminum  is  about  one  third  that  of  steel,  and  if  scantlings 
were  identical,  the  deflection  of  an  aluminum  ship,  under  a  gi- 
ven load,  would  be  three  times  those  of  a  steel  ship. 

However,  to  obtain  a  factor  of  safety  on  ultimate  strength 
equal  to  that  in  steel,  the  area  of  an  aluminum  member  in  direct 
stress,  or  the  section  modulus  of  an  aluminum  member  in  bending 
would  be  about  1.5  times  that  of  the  corresponding  steel  member 
(based  on  an  ultimate  strength  of  about  60000  psi  in  medium  steel 
and  of  about  40000  psi  across  the  welds  of  a  welded  aluminum 
structure ) . 

A  section  modulus  of  1.5  times  that  in  steel,  with  the 
same  depth  being  considered  in  both  cases,  would  result  in  a 
deflection  in  aluminum  about  twice  that  in  steel.  However,  the 
lower  modulus  of  elasticity  of  aluminum  reduces  its  ability  to 
resist   buckling  ,  and  for  the  same  factor  of  safety  both  the 
section  modulus  and  the  moment  of  inertia  in  aluminum  would  in 


b1^ 


general  be  somewhat  more  than  1.5  times  those  in  steel,  so 
that  the  deflection  of  the  aluminum  structure  would  in  general 
be  somewhat  less  than  twice  that  of  the  steel  structure. 

The  greater  deflection  would  reduce  the  natural  frequency 
of  hull  vibration,  and  sowould  affect  the  interaction  of  slamming 
stress  and  bending  moment  stress.  It  could  also  affect  shaft 
bearing  pressures.  It  has  been  proposed  as  an  arbitrary  basis 
of  design  that  aluminum  ships  be  designed  for  a  deflection  about 
1.5  times  that  of  steel  ships.  This  would  require  a  moment  of 
inertia  of  twice  the  value  of  that  in  the  steel  hull  girder, 
and  would  automatically  result  in  a  section  modulus  substantially 
more  than  1.5  times  that  in  steel,  and  a  factor  of  safety  cor- 
respondingly greater  than  in  steel. 
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SUMMARY 

Length 

= 

1000« 

Beam 

= 

al/5.75 
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B/3.3 
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= 

18000  psi 

S 

= 

10" 

Material 
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= 
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= 
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= 

18111. ll 
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= 

53387.98 

SBCRIT 

= 

48648.64 

DITOP 

= 
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= 

53. 34' 
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SUMMARY 

Length 

= 

1000' 

Beam 

= 

AL/5. 75 

Draft 

= 

B/3.3 

Depth 

= 

AL/14 

DST 

= 

18000  psi 

S 

= 

30" 

Material 

- 

ALUMINUM 

TOP 

= 

4.86" 

BOT 

= 

1.00" 

STH 

= 

18063.18 

SBH 

= 

41334.96 

SBCRIT 

= 

41021.80 

DITOP 

= 

21.72" 

DIBOT 

= 

49.71' 

STRAIGHT     CF     SMPS-ON     Thr     NEUTRAj      AXIP     LOCATION 
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SUMMARY 

Length 

= 

1000» 

Beam 

= 

AL/5.75 

Draft 

= 

B/3.3 

Depth 

= 

AL/14 

DST 

= 

18000  psi 

S 

= 

120" 

Material 

- 

ALUMINUM 

TOP 

= 

^.27" 

BOT 

= 

2.98" 

STH 

= 

17958.79 

SBH 

= 

22852.76 

SBCRIT 

= 

22768.12 

DITOP 

= 

31.43' 

DIBOT 

= 

39.99'  or  40.00' 
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Looking  at  the  following  table  we  can  also  observe  for  the 
aluminum  the  effect  of  changing  S  on  the  neutral  axis  location 
and  the  compressive  stresses  resulting  on  thebottom  plating. 


CASE 

s 

inches 

TOP 
inches 

BOT 
inches 

SBH 
psi 

SBCRIT 
psi 

DITOP 
ft 

DI  BOT 
ft 

1 

10 

5.35 

0.363 

53000 

49000 

18.09 

53.34 

2 

30 

4.86 

1.00 

41000 

41000 

21.72 

49.71 

3 

60 

4.53 

1.76 

32000 

32000 

25.83 

45.60 

*♦ 

120 

4.27 

2.98 

23000 

23000 

31.43 

40.00 

(*)  This  difference  may  result  from  inaccuracy  in  reading 
the  plot. 

Here  again  the  values  of  TOP  and  BOT  may  allow  for  a  quick 
evaluation  of  the  weight  relationship  between  the  four  cases. 


One  run  was  also  made  in  the  case  of  the  aluminum  for 
AL=500  ft  and  the  results  are  here  included. 
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top  thickness    (in.) 
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Intersection  not  accepted;  assumed  Top>  Bot  ended  up 
with  Bot > Top. 
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of» 


top  thickness  (in.) 


0.1- 


o.l- 


Qif 


O.i- 


0-4' 


o.J. 


O.0-. 


o.\-- 


o-t 


o.a     o.j     o.+    0.5- 

bottom  thickness  (in.) 
this  intersection  is  acceptable  Bot>Top 


0-6      0.7 
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SUMMARY 

Length 

= 

500' 

Beam 

= 

AL/8 

Draft 

= 

B/3 

Depth 

= 

AL/9 

DST 

= 

18000  psi 

S 

= 

30" 

Material 

- 

ALUMINUM 

TOP 

= 

0.38" 

BOT 

= 

0.6" 

STH 

= 

18089.66 

SBH 

= 

14578.14 

SBCRIT 

= 

14767.85 

DITOP 

= 

30.76" 

DIBOT 

= 

24.79" 
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HOGGING  BENDING  MOMENT=SAGGING  BENDING  MOMENT 
DESIGN  STRESS  IN  TENSION=DESIGN  STRESS  IN  COMPRESSION 

From  the  preceding  section  we  can  understand  the  need  to 
put  some  limit  on  the  value  of  the  compressive  stresses,  and  per- 
haps also  the  inefficiency  of  considering  the  value  of  S  (the 
spacing  between  longitudinal  stiffening  members)  as  an  input  to 
the  problem. 

In  this  section  we  will  consider  a  bending  moment  of  the 
same  magnitude  in  both  hogging  and  in  sagging  and  the  stress 
in  either  tension  or  compression  will  have  as  upper  limit  the 
design  stress. 

Since  the  bending  moment  is  the  same,  it  is  intuitive  that 
the  thickness  required  to  resist  tension  on  top  and  compression 
on  bottom  (buckling)  when  the  hogging  moment  is  acting  will 
also  be  respectively  the  thicknesses  to  resist  tension  on  bottom 
and  buckling  on  top  when  we  consider  the  sagging  moment.  We 
would  end  up  then  with  a  thickness  on  top  required  by  the 
hogging  moment  TOPH,  a  thickness  on  bottom  required  by  the  hog- 
ging moment,  as  well  as  a  thickness  on  top  required  by  the  sagging 
moment  TOPS  and  a  thickness  on  bottom  required  by  the  sagging 
moment  BOTS. 

Since    HM  =  SM  we  must  obtain 

TOPS  =  BOTH  =  tl  and  BOTS  =  TOPH  =  t2 

and  we  will  end  up  considering  the  largest  tl  or  t2,  to 

be  the  thickness  of  the  top  and  bottom  simultaneously. 

Assume  for  instance  that  tl >  t2  then  if  we  decide  on  tlfor 
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thickness  of  the  top  as  well  as  at  the  bottom,  we  will 

verify  that  we  have  overdesigned  because  TOPS  =  tl 

corresponds  to  BOTS  =  t2  which  plugged  into  the  expression 

for  I  gives  a  smaller  value  than  using  TOPS  =  tl  with 

BOTH  =  tl  that  will  give  a  larger  value  for  I;  this  larger 

I  value  when  used  in 

BM*Y 
STR  = 


I 
will  give  a  smaller  STR  value  than  the  design  stress  and 
we  could  use  a  smaller  thickness,  may  be  between  tl  and 
t2  for  both  the  top  and  bottom  and  satisfy  the  DST  requi- 
rement. 
This  means  that  we  must  incorporate  in  this  analysis  the  fact 
already  known  that  when  BM  (hogging)  is  equal  to  BM  (sagging) 
we  will  have  to  consider  the  neutral  axis  half  way  between  top 
and  bottom  which  means  that  the  thickness  of  the  top  shall  be 
equal  to  the  thickness  of  the  bottom  for  a  geometrically 
symmetrical  cross  section. 

Let's  determine  the  formulas  that  will  best  serve  this 
case  then: 

HM  =  SM  =  BM 

DITOP  =  DIBOT  =  D*(B*t+D*t)/(2*D*t+B*(2*t)) 

t  =  TOP  =  BOT  =  TS 


3*B*B+D*D+4*B*D 
I  =  2*D*D*t* 


D+B 

BM*DITOP    BM*DIBOT 
Now  we  can  write  STR  =  =  =  DST 
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and  the  instability  formula 

t  =  S 


*  \  0 


1-P0I-P0I 
30396* *DST 


Where  S  is  the  spacing  between  longitudinals. 

We  see  thus  that  we  cannot  have  more  than  two  equations 
and  so  we  must  seek  for  two  variables.  A  quick  inspection  and 
one  can  see  that  knowing  BM,  B,  D,  the  material  properties  and 
the  design  stress  DST  we  have  left  to  determine  what  should  be 
the  thickness  of  the  plating  (t)  and  the  corresponding  frame 
spacing  (S).  Solving  the  equations  with  this  in  mind,  we  obtain 

BM  =  1.371^3*B*T*AL*AL   lb*ft 


BM*Y 


and  from 


STR  = 


we  obtain 


BM*(D+B) 


DST 


4*D*t ( 3*B*B+D*D+^*B*D ) 


which  gives 


BM*(D+B) 


t  =• 


*J.*D*DST*  ( 3*B*B+D*D+4*B*D ) 


and  also 


S  = 


1 


0. 30396*(1-POI*POI ) -DST 
E 
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These  formulas  will  give  for  a  prescribed  set  of  dimensions 
(ft)  and  a  particular  material  subject  to  bending  BM  (lb*ft)  the 
thickness  t  (inches)  and  the  spacing  of  longitudinal  stiffening 
members  S  (inches)  required. 
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HOGGING  MOMENT  AND  SAGGING  MOMENT  OF 
DIFFERENT  MAGNITUDES 
THE  DESIGN  STRESS  IS  THE  SAME  IN  TENSION  AND 
IN  COMPRESSION 

Let's  consider  now  the  case  where  the  magnitude  of  HM  is 
different  than  the  value  for  SM. 
We  have  in  this  case: 


D*(B*BOT+D*TS) 
DITOP  = 


DIBOT  = 


2*D*TS+B*(TOP+BOT) 

D*(B*TOP=D*TS) 

SAME  DENOMINATOR 


I  =  4#D*D* 


3#B*B*TOP*BOT+ ( D*TS ) 2+2*B*D*TS* ( TOP+BOT ) 
SAME  DENOMINATOR 


HM*DITOP                         SM*DITOP 
STH  =  STS  =  


HM*DIBOT 
SBH  =  SM*DIBOT 


SBS  = 


BOT   =    \|g     \JSBH  TOP  =    \[g~  \JsTS 


where  G  =   0.30396* 


S2*(l-POI*POI) 
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Looking  at  these  equations,  we  know  that  the  result  obtained 
for  the  thicknesses  is  going  to  determine  the  value  of  DITOP, 
DIBOT  and  I.  If  furthermore  we  restrain  the  stresses  from  being 
larger  than  DST  we  will  undoubtedly  want  them  to  end  up  being 
DST  otherwise  if  they  are  smaller  we  are  not  using  the  material 
efficiently. 

We  thus  have 

DST  =  STH  =  HM*K1       DST  =  STS  =  SM*K1    where  Kl  =  DITOP/l 

DST  =  SBH  =  HM*K2      DST  =  SBS  =  SM*K2    where  K2  =  DIBOT/l 

This  tells  us  that  if  for  instance  HM >  SM 

we  should  be  considering  DST  =  STH  =  HM*K1 

and  DST  =  SBH  =  HM*K2 

because  the  other  two  equations  will  determine  lower  stresses 

and  are  thus  within  the  requirement  of  not  exceeding  the 

design  stress,  being  automatically  satisfied. 
Now  from  these  equations  we  can  see  that  if  the  design  stress 
required  for  the  top,  (in  tension  for  the  case  of  HM>SM;  note 
that  for  the  case  SM  > HM  the  other  two  equations  should  be  con- 
sidered and  the  conclusions  are  similar),  is  going  to  be  the 
same  one  as  required  for  the  bottom  (  in  compression  ),  we  must 
have  Kl  =  K2  and  so  DITOP  =  DIBOT  and  the  neutral  axis  will 
again  end  up  being  located  halfway  between  top  and  bottom  for 
a  geometrically  symmetrical  cross  section. 

We  can  here  use  the  same  two  equations  derived  in  the 
previous  section  to  determine  (t)  and  (S).  The  moment  to  be 
considered  here  will  be  the  larger  of  HM  or  SM  instead  of  BM. 
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HOGGING  MOMENT  LARGEH  THAN  SAGGING  MOMENT 
DESIGN  STRESS  IN  COMPRESSION  LARGER  THAN  DESIGN  STRESS  IN  TENSION 


Thesituation  we  are  analysing  in  this  section 

HM>SM 
DSC  >  DST 
is  perhaps  more  interesting  than  the  previous  ones  because  its 
study  provided  us  with  a  much  simpler  method  for  determining  the 
thicknesses  and  the  neutral  axis  location;  we  call  here  the  atten- 
tion of  the  reader  to  fully  understand  the  steps  of  this  section 
because  the  following  sections  till  the  end  of  this  work  will  use 
this  same  method  of  analysis  applied  totheir  own  particular  case. 

Indicated  in  the  next  diagram  are  the  ranges  within  which  we 
may  stress  the  material  in  order  to  satisfy  the  requirement  of 
not  exceeding  DSC  or  DST. 


— DSC    *■ 

p* DST-*- 

i 

:> 

COMPRESSION 

TENSION 

Let's  consider  now  that  we  will  try  to  determine  the  best 
position  for  the  neutral  axis  by  first  gaining  some  insight  into 
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the  geometry  of  the  stress  distributions  and  the  stress  limits 


already  defined, 


DSC- 


./ 


DITOP 


DIBOT 


Acceptable  SBH 


DST- 


Unacceptable 
too  large  STH 


/ 


NEUTRAL  AXIS 


I 


From  the  previous  diagram  we  may  easily  understand  that 
if  we  decide  to  stress  the  material  to  its  maximum  in  tension 
and  compression  when  under  the  effect  of  a  sagging  bending 
moment  (SBS  -  STS)  we  end  up  with  a  neutral  axis  closer  to  the 
bottom  due  to  the  fact  that  DST <  DSC  and  we  see  that  when  consi- 
dering the  hogging  moment  the  only  possible  alternative  left 
is  to  limit  the  stress  caused  by  this  moment  in  tension  on  top 
STH=DST;  in  this  case  the  compressive  stress  on  the  bottom  SBH 
would  automatically  fall  inside  the  acceptable  range. 

If  for  the  hogging  moment  we  choose  to  allow  the  stress 
in  compression  on  the  bottom  to  attain  its  maximum  value  SBH=DSC 
we  see  from  thediagram  that  we  are  forcing  the  STH  to  fall  out- 
side of  the  acceptable  range.  We  are  left  in  this  situation  with 
an  acceptable  solution  so  far  and  we  must  investigate  it  further. 
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At  this  stage  we  have 


I  qT.q=nsn 

STH= 

DST 

1                           \k 

I        \ 

/ 

/ 

/ 

/ 

N.A. 

I         / 

\     I 

sbh/ 

SBS= 

^DST 

We  may  see  from  the  diagram  that  we  are  asking  for  SBS=STH- 
=DST  to  have  the  same  value  and  at  the  same  time  to  have  the 
neutral  axis  located  closer  to  the  bottom  DITOP>  DIBOT. 

We  note  that  the  equations  that  give  STH  and  SBS  are 


and 


STH  =  HM*DITOP/l 
SBS  =  SM*DIB0T/I 


and  since  I  is  the  same  STH  can  never  be  of  the  same  magnitude 
as  SBS  because 

STH  =  HM*DITOP/l  >  SBS  =  SM*DIBOT/l 


both  because  DITOP  is  being  set  larger  than  DIBOT  and  also 
because  we  are  considering  HM  > SM  to  begin  with.  Thus  we  are 
left  with  the  alternative  of  making  DIBOT  larger  than  DITOP 
to  make  it  possible  to  satisfy 

STH  =  HM*DIT0P/I  =  SBS=  SM*DIBOT/l 
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Thesituation  would  then  be 


DSC 

vSTS 

STHf 

=  DST 



— 

4 

/    i 

/ i 

... 

\        i 
\      i 

\  j 

H..A,. 

DSC= 

5BH 

SBS^ 

--DST 

Besides  this  is  only  natural  because  what  we  are  doing  here 
is  to  make  the  material  to  be  stressed  to  its  allowable  limit 
while  resisting  the  larger  of  the  two  moments;  hogging  in  this 
case.  This  situation  will  be  worked  out  here  as  an  example  and 
it  seems  to  be  quite  clear  that  similarly,  cases  where  SM>  HM  or 
DST>DSC  or  other  situation  like  DST>DSC  with  HM  >  SM  or 
DST<DSC  with  SM  >  HM  could  be  solved  applying  a  similar  reasoning, 
So,  having  HM>SM  and  DSC >  DST  we  should  consider  the  smaller 
design  stress, in  this  case  DST,and  have  the  material  stressed 
to  this  upper  limit.  Why?, because  if  the  larger  is  considered, 
DSC  in  this  case,  there  is  no  way  we  can  have  the  material  stre- 
ssed to  this  higher  limit  in  both  hogging  and  sagging  while  con- 
sidering the  smaller  we  may  attempt  to  do  so.  This  can  be  better 
understood  from  the  two  next  diagrams. 
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DSCtSTS 


DSTrSTH 


DSCfSBH 


oo  large 

SBS 


DSC=STS 

DST, 

Tc 

1    \ 

!      .    \ 
1              \ 

!                  \ 

i                X 

i 

1 

1  y 

/S    i 
s          i 

N.A. 

1                / 

1                       yS 

\           / 

X      ! 

\i 

DSC=SBH 

DST=t5BS 

Toolarge 
TH 


So,  considering  the  lower  design  stress,  in  this  case  DST, 
we  will  stress  the  material  to  this  maximum  in  both  hogging  and 
sagging  thus  giving 

STH  =  HM*DIT0P/I  =  DST 
SBS  =  SM*DIBOT/I  =  DST 
Note  here  that  we  already  know  that  the  neutral  axis  will. 
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be  closer  to  the  top  DITOP<DIBOT  which  means  that  we  must  obtain 
a  thickness  for  the  plating  at  the  top  larger  than  at  the  bottom  and. 
so  TS=TOP  which  means  that  we  will  be  using  the  following  for- 
mulas 

D*(B*BOT+D*TOP) 


DITOP  = 


P 
D*TOP*(B+D) 


DIBOT 


P 
4*D2*T0P*( ( 3*B2+2*B*D)*BOT+(D2+2*B*D)*TOP) 


I  = 


where   P  = ( 2*D+B)*T0P+B*B0T 
we  have  also 

STH  =  DST  =  HM*DIT0P/I 

SBS  =  DST  =  SM*DIB0T/I 

SBH  =  DSC  =  HM*DIB0T/I 
From  these  three  equations  we  may  combine  the  two  first  to  obtain 

HM*DITOP  =  SM*DIBOT 
and  combining  the  second  and  third  we  obtain 

DST*HM  =  DSC*SM 
After  combining  these  two  equations  we  obtain 

DST*DIBOT  =  DSC*DITOP 

This  equation  could  have  been  obtained  more  easily  just 
by  looking  at  the  diagram  we  are  considering  and  notice  the 
similar  triangles 


1^2 


DSC 


Note  also  that  what  we  obtain  from  this  equation  is  to 
determine  the  N.A.  location  by  making  DITOP  and  DIBOT  proportio- 
nal to  DST  And  DSC. 

Proceeding  with  this  equation  we  obtain: 

DST*(D*TOP*(B+D))/P  +  DSC*(D*B*B0T+D2*T0P)/P 

TOP*(DST*B+DST*D-DSC*D)  =  DSC*B*BOT 

and  making 

AC  =  DST/DSC+D*(DST-DSC)/(DSC*B) 

we  obtain 

BOT  =  AC*TOP 


We  need  now  an  equation  to  stress  the  material  to  its  design 
limit  and  this  can  be  done  by  either  using 

STH  =  DST  =  HM*DIT0P/I 
or    SBH  =  DSC  =  HM*DIBOT/l 

Prom  these  two  equations  I  prefer  to  use  the  second  especially 
because  it  will  simplify  the  calculations  that  follow  as  can 
be  understood  from  the  simpler  expression  that  gives  DIBOT  as 
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compared  with  the  one  that  gives  DITOP. 

We  will  then  consider  here 

SBH  =  DSC  =  HM*DIBOT/l 

which  gives 

HM*D*TOP*(B+D) 

DSC  = 


4*D2*T0P( (3*B2+2*B*D)*BOT+(D2+2*B*D)*TOP) 


Simplifying  we  obtain 

»HM*(B+D)             3*B2+2*B*D 
TOP  = *BOT 

DSC*4*D* ( D2+2*B*D )    D^+2*B*D 


or  making 
A  = 
AB  =  (3*B2+2*B*D)/(D2+2*B*D) 


A  =  HM*(B+D)/(DSC*4*D*(D2+2*B*D)) 


we  have 

TOP  =  A-AB*BOT 
replacing  here 

BOT  =  AC*TOP 

we  obtain 

A 
TOP  = 


1+AB*AC 


Having  already  found  TOP  and  BOT  we  may  use  BryanS formula  to 
determine  S; since  the  thickness  of  plating  is  for  this  case 
smaller  at  the  bottom  we  must  consider  this  thickness  and 
along  with  it  DSC  because  this  will  be  the  larger  compressive 
stress  the  bottom  will  be  subjected  to. 


lkk 


E 


BOT 


From  DSC  =■ 


0.30396*(1-POI*POI) 


We  obtain 
S  =• 


\ 


•*BOT 


1-POI--POI 

0.  30396*— *DSC 

E 


or  making 
Z  =■ 


v 


1-P0I*P0I 

0.  30396* —-"DSC 

E 


we  have 

S  =  Z*B0T 

A  computer  program  was  written  and  applied  to  two  parti- 
cular sets  of  dimensions  and  two  different  materials  namely  steel 
and  aluminum.  The  results  are  presented  next  in  the  computer 
printouts. 
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HOGGING  MOMENT  AND  SAGGING  MOMENT  OF  EQUAL  MAGNITUDE 
DESIGN  STRESS  IN  TENSION  AND  DESIGN  STRESS  IN  COMPRESSION  OP 
EQUAL  VALUE;  LATERAL  LOAD  ALSO  CONSIDERED 

The  lateral  load  being  here  considered  will  be  a  "head" 
of  salt  water  of  eight  feet  on  the  top  plating  (HTOP=8  ft) 
and  a  "head"  of  salt  water  equal  to  the  draft  on  the  bottom 
plating  (HBOT=DRA). 

We  will  use  here  the  simplified  approach  as  suggested 
in  the  following  note;  this  simplified  approach  is  at  least 
conservative . 

Note:  Bleich,  from  his  experiments  and  theoretical  studies 
on  the  effect  of  normal  pressure  in  the  buckling  strength  of 
thin  panels,  concluded  that  the  buckling  strength  is  increased 
when  combined  with  a  sufficiently  large  normal  load,  and  that 
the  increase  is  greater  for  thesimply  supported  edges  situation 
than  with  fixed  edges.  For  heavier  plating  ,  where  thedef lection 
does  not  exceed  half  of  the  plating  thickness,  the  increase  in 
buckling  strength  resulting  from  normal  load  is  negligible. 
It  was  therefore  suggested,  for  simplicity,  that  the  critical 
stress  be  computed  for  in  plane  loading  only. 

The  tertiary  stresses  we  are  interested  in  considering 
here  are  those  exerted  at  points  B  directed  longitudinally;  these 
will  combine  with  the  longitudinal  bending  stresses  or  primary 
stresses  as  described  in  a  previous  section. 

Since  now  we  have  DST=DSC  and  HM=SM=BM  we  may  think  of 
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a  diagram  for  the  stress  distributions  as  indicated  next  figure: 
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SBS          ■,ST33B 

COMPRESSION    TENSION 


From  this  diagram  we  may  understand  that  the  neutral  axis 
should  be  located  closer  to  the  bottom  which  means  that  the  thi- 
ckness at  the  bottom  will  be  larger  than  at  the  top  BOT>TOP  . 
So,  we  may  already  write  the  following  formulas  where  TS  was 
replaced  by  BOT 

DITOP  =  D*(B+D)*BOT/P 

DIBOT  =  D*(B*TOP+D*BOT)/P 

I  =  ij.*D*D*(  (3*B2+2*B*D)*T0P+(D2+2*B*D)*B0T)*B0T/P 


where 


P  =  2*D*B0T+B*(T0P+B0T) 
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We  may  now  think  of  the  following  relationships  to  deter- 
mine TOP  and  BOT. 


DITOP 


N.A. 


DIBOT 


SBH 
STH*DIBOT  =  SBH*DITOP 

or  knowing  that 

S   2 
STR3T  =  0.152222*HT0P*( )* 

TOP 

S  9 
STH3B  =  0.152222*HB0T*( )£ 

BOT 


and  anticipating 

S  =  Z*TOP 
where         


\ 


1-P0I*P0I 
Z  =  l/(  \  0.30396* *DST 


we  have 

2     S  9 
Z2  =  ( )2 

TOP 

S  2    9 
If  we  now  make  ( )  =  Zc 

BOT 
in  order  to  simplify  the  equations,  we  see  that  when  we  do  this, 
since 
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->2<  (— )2 


BOT 


TOP 


we  are  increasing  the  value  of  STR3B  and.  this  means  that  since 

DST  =  SBH+STR3B 
we  are  decreasing" slightly"  the  value  of  SBH  that  could  be  allo- 
wed to  reach  a  larger  value;  so  ,  the  substitution 


( )2  =  z2 

BOT 


is  on  the  safe  side  and  has  the  advantage  of  simplifying  the 
equations  considerably. 

So  we  may  write 

STH  =  DST-STR3T 

SBH  =  DST-STR3B 

STR3T  =  0.152222*HTOP*Z2 

STR3B  =  0.152222*HBOT*Z2 

A  =  DST-STR3T 

AB  =  DST-STR3B 
so 

A*DIBOT  =  AB*DITOP 
and  we  obtain  replacing  DITOP  and  DIBOT  by  their  equations, 
and  simplifying 


AB*(B+D)-A*D 

TOP  = *BOT 

A*B 


or  making 
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AC  =  (AB*(B+D)-A*D)/(A*B) 
we  finally  obtain 
TOP  =  AC*BOT 


Now  in  order  to  stress  the  material  to  its  design  stress 
we  may  either  consider  STH,  STS,  SBH,  SBS. 

Due  to  ease  of  calculation  we  may  either  choose  STS  or 
STH;  the  expressions  will  be  exactly  the  same.  I  will  consider 
here 

STH  =  BM*DIT0P/I 
and 

STH  =  DST-STR3T  =  A 
Substituting  for  DITOP  and  I  we  obtain 


A  =• 


BM*D*(B+D)*BOT 


i+*D2*B0T*  (  (  3*B2+2*B*D )  *TOP+  ( D2+2*B*D )  *BOT ) 


and  after  simplifying  we  arrive  at 


BM*(B+D) 


BOT  = 


3*B2+2*B*D 


A*4*D*(D  +2*B*D)    D2+2*B*D 


■if- 


TOP 


Making 

AD  =  BM*(B+D)/(A*^*D*(D2+2*B*D) ) 
AE  =  (3*B2+2*B*D)/(D2+2*B*D) 

we  obtain 

BOT  =  AD-AE*TOP 
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substituting 

TOP  =  AC*BOT 
we  obtain 


AD 


BOT  =• 


1+AC*AE 


Considering  now  the  instability  formula  (BRYAN)  we  notice 
that  since  the  top  will  have  the  thinner  plating  we  will  consi- 
der TOP  stressed  to  the  maximum  allowed  DST. 
So  we  start  with 


Z  =  l/( 


\l 


1-P0I*P0I 

0.30396*- *DST  ) 

E 


S  =  Z*TOP 

The  total  stresses  will  then  be  given  by 
TSTH  =  BM*DIT0P/I+STR3T 
TSTS  =  BM*DIT0P/I+STR3T 

TSBH  =  BM*DIBOT/l+0.l52222*HBOT*S2/BOT2 
TSBS  =  BM*DIB0T/I+0.152222*HB0T*S2/B0T2 


The  computer  printouts  for  the  same  cases  analysed  in  theprevious 
section  are  presented  next. 
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strfnght  of  ships-on  thf  neutral  axis  location 
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STRENGhT    OF    SHIPS-ON    THE    NEUTRAL    AxIS    LOCATION 
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Note  here  that  because  we  have  in  this  run  BOT>>  TOP  the 

approximation  made  to  simplify  the  formulas,  namely 
S 

Z  = 

BOT 

becomes  very  inaccurate  compared,  to  what  the  actual  value 

should  be 

S 
Z  = 

TOP 

this  is  the  cause  for  obtaining  such  low  stresses  at  the  bottom. 
We  might  refine  the  result  for  this  case  as  is  done  for  the 
"1000  ft  steel  case"  also  included  in  this  section  which  presen- 
ted still  another  difficulty.  Since  all  steps  to  have  the  bottom 
stressed  to  DST  are  explained  in  detail  for  such  case  they  will 
be  omited  here. 
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This  case  gave  a  result  close  to  the  correct  solution  with 
the  bottom  stressed  to  almost  DST. 

Next  is  presented  the  "1000  ft  steel  case"  where  still 
another  difficulty  was  encountered  and  a  procedure  to  have  the 
bottom  stressed  to  DST  is  exemplified. 


STRENGmT   of   ships-qn   the   neutral   axis  location 
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(,  VIC  OPERATING  SYSTEM  VERSION  1  REVISION   012  P3/0W74  GENERATED   04/30/74 


3B 
?5 
>t 

[*0 

il 

^0 


HANDLING  CHARGE 
LINES  PRL\TED  FR2 
CARDS  RE  A;; 
PLOTTER  VECTORS 
MODEL  73  SECONDS 
vCDEu  80  SECONDS 


,3b 
1.25 
1,50 

.25 


*  2  5 . 0  0 
512.50 


JOB 
K  LM 
K  CD 
1  7  0  7 
HOUR 

wruR 


TOTAL  CHARGE  * 


.35 
•16 

.13 

,cv 

.  15 

.00 

.79 


KIR  hS?     1^731   LOGGED  OUT  ?4/3?/7^ 


27.^7  LEFT  AFTER 


5J  LOGI  JS. 


174 


It  does  not  come  as  a  surprise  that  the  program  does  not 
work  for  the  "1000  ft  steel  case"  since  in  this  situation  we 
obtain  STR3B  DST  which  cannot  be  accepted  and  also  this 
program  is  not  prepared  to  handle. 
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U  =   W  =  15.182*10"3 

z  =  i/u  =  65.867 

HB0T  =  DRA  =52.7 

STR3B  =  0.152222*HB0T*Z2  =  34803.57    DST 
Since  we  understand  that  the  program  should  work  for  DST>  STR3B 
we  make  the  next  run  with  DST  =  35000  to  prove  that  in  such  case 
a  solution  is  obtained^and  after  proceed  to  obtain  the  actual 
solution  for  the  case  we  were  considering  with  DST  =  25000. 
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Now  let's  consider  our  situation  with  DST  =  25000.  Once 

the  situation  is  such  that  due  to  the  approximation  made  in 
S 

=  Z 

BOT 
we  end  up  with  a  STR3B>DST  we  must  consider  a  value  Zl  lower 

than  Z  to  bring  STR3B  down.  This  can  be  also  understood  since 

it  is  acceptable  to  make  here 
S 

=  Z 

TOP 
but  because  we  obtain  with  the  program  as  is  T0P«B0T,  making 

the  approximation 
S 

=  Z 

BOT 
makes  Z  attain  a  much  larger  value  than  it  actually  should. 

We  must  then  consider  Zl<  Z  =  65.8  and  may  start  by  trying 

Zl  =  35 

The  computer  printout  in  this  case  is  as  follows. 
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This  time  we  obtained  a  realistic  but  still  unacceptable 

result;  we  notice  that  the  bottom  is  being  stressed  to  31000  psi 

when  it  can  only  go  up  to  25000  psi;  besides  we  assumed 
S 

Zl  =  . =  35 

BOT 
and  end  up  with 

S    225.^ 

= =  44 .  46 

BOT    5.07 
This  means  that  the  Zl  value  should  be  higher. 

We  may  now  try  Zl  =  4-0;  the  computer  results  for  this  case  are 

shown  next. 
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In  this  case  we  assumed  Zl  =  40  and  ended  up  with 

S    199.5 

BOT 


=  15.09 


6.59 

and  a  stress  at  the  bottom  of  only  19000  psi  which  is  too  small, 

We  may  now  try  to  use  a  better  Zl  value  and  can  do  the 

following: 
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we  will  use  this  Zl  value 
^   The  results  obtained  with  this  Zl  value  are  the  following: 
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The  results  show  that  we  have  obtained 


213.9701 


Zl  =• 


BOT    5.7^18 


-=  37.265 


while  it  was  assumed 

Zl  =  37.617 
We  could  further  correct  the  result  until  Zl  assumed  and  Zl  actual 
would  coincide  but  it  seems  that  there  is  no  need  to  refine  more 
the  solution  for  the  purpose  here  involved. 


Note  that  this  procedure  should  also  be  followed  for  the  other 
cases  in  order  to  bring  the  stress  at  the  bottom  to  its  maximum 
allowed  DST  for  such  cases  as  well. 
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This  1000  ft  aluminum  case  shows  again  a  stress  at  the 
bottom  much  lower  than  DST.  It  should,  be  corrected  as  done  for 
"1000  ft  steel"  in  the  previous  situation. 


195 


HOGGING  MOMENT  LARGER  THAN  SAGGING  MOMENT 
DESIGN  STRESS  IN  COMPRESSION  LARGER  THAN  DESIGN  STRESS  IN  TENSION 
LATERAL  LOAD  ALSO  CONSIDERED 


With  the  help  of  what  was  understood  for  the  similar  si- 
tuation without  the  lateral  load  being  included,  we  may  here 
proceed  at  once  to  look  at  what  the  diagram  of  stress  distribu- 
tions might  look  like.  Actually,  the  diagram  that  will  be  con- 
sidered first  is  basically  the  same  with  only  a  small  change 
due  to  the  addition  of  the  tertiary  stresses;  since  these  will 
probably  be  of  the  same  magnitude  on  both  top  or  bottom  and 
assuming  that  their  actual  value  does  not  differ  considerably, 
what  we  may  start  looking  into  is  a  diagram  as  follows: 


DSC 


DST 


We  would  consider  here  the  same  basic  distribution  of  longitu- 
dinal bending  stresses  thus  obtaining 
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We  will  then  expect  to  obtain  TOP>  BOT  and  DITOF<  DIBOT 
and  we  may  already  consider  the  respective  expressions  for 
DITOP,  DIBOT  and  I  with  TS  replaced  by  TOP 

Prom  the  diagram  we  may  consider  thefollowing  equations 

STH*DIBOT=SBH*DITOP 
and  notice  here  that  if  we  want  to  stress  the  material  to  the 
allowable  design  stresses  we  have 

STH  =  DST-STR3T 

SBH  =  DSC-STR3B 

we  may  combine  these  equations  to  obtain 

(DST-STR3T)*DIB0T  =  (DSC-STR3B)*DIT0P 

STR3T  and  STR3B  are  given  by 

S   2 

STR3T  =  0.152222*HTOP*( -) 

TOP 


STR3B  =  0.152222::-HB0T*( \ 

BOT 
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anticipating  here  that  S  =  Z*BOT  because  BOT  is  the  thinner 
plating  we  have 


S  .  o    o 

( — r=  zd 


BOT 


and  thus  may  write 

STR3B  =  0.152222*HB0T*Z2 


Since 


( )2<(— )2 


TOP     BOT 
because  TOP >  BOT,  if  we  make 

( )2=  Z2 

TOP 

to  simplify  the  equations  that  will  follow  we  will  understand 
that  in  doing  so  we  are  only  increasing  the  value  of  STR3T. 
This  as  we  see  fron 

STH  =  DST-STR3T 
will  make  us  have  STH  with  a  slightly  lower  value  than  the  one 
that  is  allowed,  so  we  are  on  the  safe  side  when  we  consider 

( )2  =  Zd 

TOP 

Now  if  we  end  up  too  much  on  the  safe  side  and  we  have  a  TSTH 
too  low  we  may  bring  it  to  a  value  closer  to  DST  by  considering 
a  Zl  value  as  done  and  exemplified  in  the  previous  section  and 
done  again  in  this  section  for  the  case  1000  ft  steel  where 
for  this  particular  exemple  it  is  TSBH  that  needs  correction 
as  explained  in  detail  . 
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We  may  then  write 

1 


Z  = 


\ 


1-P0I*P0I 

0.  30396* *DSC 

E 


STR3T   =    0.152222"HT0P*Z2 

STR3B  =  0.152222*HB0T*Z2 
and  from 

(DST-STR3T)*DIB0T  =  (DSC-STR3B)*DIT0P 
we  have  if  we  make 

A  =  DST-STR3T 

AB  =  DSC-STR3B 

P  =  (2*D+B)*T0P+B*B0T 

D*T0P*(B+D)       D*(B*B0T+D*T0P) 

A* — — —  =  AB*- _______ 

P  P 

which  gives  after  simplification 

A*((B+D)-AB*D)*TOP  =  AB*B*B0T 
and  making 

AC  =  (A*(B+D)-AB*D)/(AB*B) 
we  have 

BOT  =  AC*T0P 

Following  th  e  suggestions  of  a  previous  section  we  will 
also  use  the  equation  that  will  make  the  bottom  to  be  stressed 
the  most  in  compression  under  the  larger  of  the  moments  being 
considered  (HM  in  this  case). 
We  obtain 

SBH  =  HM*DIB0T/I  =  AB 
Replacing  DIBOT  and  I  by  their  expressions  we  obtain 
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HM*D*TOP*(B+D) 

AB  = 


4.*D2*TOP*  (  (  3*B2+2*B*D )  *BOT+  (D2+2*B*D )  *TOP ) 


simplifying  we  end  up  with 

HM*(B+D)         3*B2+2*B*D 

TOP  = k *BOT 

AB*4-*D*(D2+2*B*D)    D  +2*B*D 


now  making 

AD  =  HM*(B+D)/(ABJ<^*D*(D2+2*B*D)  ) 

AE  =  (3*B2+2*B*D)/(D2+2*B*D) 
we  have 

TOP  =  AD-AE*BOT 
substituting  here 

BOT  =  AC*TOP 
we  obtain  finally 

AD 

TOP  = 


1+AC*AE 


Considering  here  again  the  bottom  as  the  location  of  the  thinner 
plating  (in  this  particular  situation  we  are  looking  at),  and 
knowing  that  it  will  be  subject  to  a  compressive  stress  of 
magnitude  as  large  as  DSC  we  may  again  write   S  =  Z*BOT 
We  finally  obtain  the  total  stress  on  top  and  bottom  in  both 
tension  and  compression  by  adding  the  longitudinal  bending 
stresses  to  the  tertiary  stresses. 

TSTH  =  HM*DITOP/I+0.152222*HTOP*S2/TOP2 

TSTS   =    SM*DITOP/I+0.152222*HTOP*S2/TOP2 

TSBH  =  HM*DIB0T/I+STR3B 

TSBS  =  SM*DIB0T/I+STS3B 
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NOTE:  since  usually  the  secondary  stresses  are  assigned 
a  constant  value  like  2000  psi  or  3000  psi  they  can  be  a- 
ccounted  for  when  giving  a  value  to  the  design  stresses. 


201 


STRfnOht  OF  SWIPS-ON  TWf  *FUTRA|   •'-  X I  c-    LOCATl 


I  i 


I 


c 

c 

c 

c 

c 
c 
c 
c 


c 
c 


*AL*  IS  THF.  L^h  •  .t 
\L  =  :••,/:». 

i  =  A  L  /  8 . 

,  OR  A  *  I  '•   T  t ,  K  ->  '„■  :  P  F 

r)RA  =  B/3« 

i  =AL/9« 

*POI*  IS  rnK  o\    I      '  s  katio  OF  IHi  MATERIAL 

,-.  >  ALU  f r<  '•'■" 

■'     I  =.  •  3-- 

*i  »     I?    IN!     Yr  iNi    i-    MOnUL  IS    OF     rHE       &TERIAL 

FOR     ALU     iM'" 

-  =  1  ?  .  »  l     .  »  -  < 


.     -7  •     T        T|i    ,iE!   !Gi  !     STKt  3S     H 


<     -cib 


TON 


=5  T  =  ?  5  fi    i    . 

*r)sc*  is  the  -5F<  rG'  Stress  if  rOMp  ession 

>C=3      '      '  • 
-H"  *      IS     THf-       :-i';,  ■  I    I!  rNUlNG     Mr)     ENT 

|M=1  .37  143*1?  *^~<      A  | 

»o:-»     IS     THJ       i  ■  ,v    ,  F NO INC     Mn     cl  J 

r  f  =  R  ¥QR  A  *  AL*'\| 

.  =.    . '-     3°6*(1*«P!    [*P^1)*USC/E 
i   = c      -  1  (      ) 

Z= 1 • / u 

jrHTOP*     IS     THE     hFaD     OF     *ATEF     ON     UFO 

hT(  P  =  8« 

.r::rT»     IS    THE'     M F. A d       r     WATER    0 '■      F  0  T T 0 M 

Hri01  =!  ^ 

S  r  -  3  1  =•  ■  I5?2?a*t-T  V  *l,.7. 

STR3W  =  /.ibp2??*--.lOT*Z^Z 

A=fST-STR3T 

AH  =  CSC-5iTF3H 

A  C  =  (  A  *  (  •    f  f  ■  )  -  A  3  *  '     )/(A?,,RI 

> .J*  =  h    +  (     +c,)/(^n^^««ri*«f:. » r)  +  p .  *  r  *  o )  ) 

\  t_  =  (  ?  .  *      *  r-  +  2  .   v  -    -       |  /  (  i     .,+■>.»  I  -  *  •-   ) 

*rOP»     I'      THE     rMlC<'  ='':    UF     TM|      rUP    t-LATlNG 

«-  CT*     I  :•     THE     rHlC<vpS^    OF    THE     liOTT     1    PLATING 

tUF =AD/ I  1 • +AC « A  .  ) 

hOT  =  AO  rop 

.'•  TTt  c?m:  )  rO's 
! ?     pORMA T (  '  I  •  ,  ^     / »  ■  rOP=  » c IP • M  ) 

.-'- ITf  (  ^,  2V  )  HOT 
?C    p  )RMAT  (  '     •/?/<,«  r,Tr=  'c-1     .  <t  I 

*S«     IS    THE     SP^CI     "5    nf     LOMGITUTI'JAL    STIFFENING    f-EMBERS 

S  =  Z  *  B  0  T 

■  - 1  t  e  (  5  >  ? i  i  s 

T--      pORMAT(   '       •   1  £  3  /  ,   1    ;3  1  n    ,  .  ..   ) 

p»p «  *0*T  OP  +  b  * ( T     •  +  - 1 T  1 

.    IIT0P*     rS     Trlr     "MST        -(■     )p     thf     T(  p    FRO"    THP     NELTRAL     A  <  T  5.' 

)ITOP=D» \     •        T+"i*T"DD)/p 

••  Rite  (  ^»  a;  )f  1  • 

/a  ••     FORMAT  (  '      '  /  1      v  ,  •    »I  Tr)F  -  «  F  J  ?  .  <t  ) 

*CIBOT*     IS     THc        '  "  r  ■.      t       IF     THF    BOTTOM    FROM     THF    SF.\T     ,\l_     *'■  .  I   ; 


202 


STRFNtVT        ^      .  i- J  PS-Ot<     T  ,- 


j  r 


'.  <  T1-    LOCAl  1 '  N 


• 


c 
c 
c 
c 
c 
c 


1 
11 

z 

A 

s 
Dl«rT 

STCPiT 
TSBS 

•  v 


If  ' 

•  IT 

4        .     -     I 
.  I* 

-  rsi 

.  i  -  r 

«  T  Si- 

tsts 

- 

. '« n 
•' 

.   -IT 

.-.  ':  I  T 

F  "    i  ■ 

T  ■  • 

■  H  fT 
P  0  K » 

"  IT 
r  "  -  ■• 

FOov 
f.  ■  N  0 

■  66  - 

,.  n    • 

f,    .  (, 

<  c: 
•'  -> 

■7U 

?  2  q  . 
.  37 
<t| 

74C 


T=D*< 

I  (  5  •  • 
&  T  (  • 

.        T 

if  •  *  I     ■• 

H  *  I S 

c,  „  T  '. 

• 

5  *  I S 

-  I"  • 

I 
=  -<■•-  r 

=  <3l I '  r 

\  r=TC 
!  (5/6 
\T(  ' 

:  i  -;  -  ; 

ft  T  |  ' 
E  (  5  }  - 
ft  T  {  ' 

ssHV  - 
=  <-'••   - 
IT=     " 
E.  (  b '  S 
AT(  ' 
E  (  5  i  1 
ft  T  <  ■ 
E  {  5  /  1 
ft  T  (  ' 

CV] 
CV] 

C  v  ] 
CV] 
C  \l  ] 
[V] 

CLD 

CL] 

Ci_j 
[  -] 


ll 
';  1 

C  .  ( 

ri 

T.., 
Tr-f. 
T  ..- 
IS 

:c^ 
I j\ 

ITU 

P.  l 

i  I 

'  ,  1 

)  ' 

•  .  ' 

7   )  S 

'  ,  3 
[ 

I     ' 
T»H 
v    )  1 
1  .  1 

;"■  j 

',! 

1 1: ) 

•  ,  i 


1  u 

- 

■ 

^7 
^T 
!  i 

•  - 
I    It: 

•   /  i 


» T  a  P  I  /  *- 


I 


;~    ( 

:  1     ' 


■     > 

•  + 
/i 


( ri 

I  -i 


-    i 

*     ■ 
i 

I      T 
:      '  i 

•AL 
•AL 

;22 

.  22 


v-  • 

1     I 

'I.     ) 

IF 

n 

31 


2  2  ■ 
6  * ! 


f- 


'+  ) 
TlA 

*jn 
Gup 

51   E 

Dl  " 
)l  £ 
■  F  •; 
<■■'- 

u  - 

h  t  n 

) 


>  'f-    i  HE    v    sr;r  noi     j 
o*Po-r  +  (n*o+2.  "R * r  )  * t of 

r(  HOGGING 

rr  ftGGING 

Tf  |  QGGING 

Tr  - AGGING 

f  \    rop 

53    C   v     rtOT 

n«;:»'-/(T  n'*Tr  pi 

;*S*<V(TOP*TCP) 
(t«-FUl*POl)) 


-  /  . 
P,  /: 


''    . 

nT* 

T,.     , 

f  >, 


1  '  5  r-..=  i  F  i  .• .  3  ) 
■  ■  I       i  r  =  '  f  1 1  •  3  ) 

-  +  '.-  t  -  ->  b 

/  t     .  3  ■»  3  S  6*F*S*(lt«iPai*P0I) 

"  i  s  ^  i  =  ,  r  1.  j « 3 ) 

*  T 
•  ■:  '  •  :  ■  I  r  -  '  p  i  .    .  3  ) 


iL 

S  ;  - 

Cwj 

a 

(  67  :  CV] 

DR< 

•> '.  3      [  V  - 

n 

c 

['J 

.R 

I  6'JC  [V] 

DST 

■  «b  A  h  C  V  - 

DSC 

*■*  ^<* 

6  js 

C-0 

•i 

?6C»  r  v1 

U 

-^      [S: 

'      M 

HToF 

.  ^  -T 

.  •  /  J 

HBOT 

0  6  r  '♦  r  v  ] 

STf  3T 

F56DC  CV; 

A  B 

■^•-;. 

Cw3 

AC 

PfeFCfV] 

AP 

?,  7 1      C  V  ] 

/  F 

■J  C  7 

p,>- 

C|  ] 

10 

P0C-'H  CS] 

(■    I 

■  2  5  --  r. L ; 

c 

-;  - 

7?" 

Cv3 

p 

0  7  F  4  r  V  ] 

DITOP 

".  3 1  <».  C L  ] 

/. 

3 

7^r 

[  >'  j 

3MI 

0  7 ;  •  3  C  V  ] 

7  5iT  H 

?-73CCV; 

TSTS 

6 ', 

•  5  •  s 

''  3 

7t' 

053c  CL] 

80 

■74 ;,  TV; 

t  ;-"-■ 

?,H-f-;  I- 

gnr 

'"'.j 

90 

0  6 1  2  C  L  ] 

10f 

1  h  **  ■■'  C  L  1 

1  1 

OGPAM    L*rFl-^: 

?K70      »  v  A  I  'I  - 

2    56 

• 

P  r  ^  c, 

2D4E     •7c"f:T 

2C 

• 

■ 

?BB6    Aj>;T 

28L2 

3     ?T 

2D52    .n 

2E?8 

1 0 

IFi?  ' 

try-poimts: 

?7C0    .p 

28    ■ 

.  a 

2332 

•  r  p 


15  ARG 


p<  AV 

/  l  n  g 

2E6<» 

si 

P7CH 

• 

2r  a/ 

.  '  1 E  3 

2C2F 

.  i 

2  A 

L 

57  4  A 

.5 

2D?C 

■  '■ 

2  y  S  ;j 

t 

'2'-,AK     /Lr'~i 


2AJH     f:1*^ 


/iujfa     A  I 


203 


rop=  '. t  -    •< 

DT=  .'.->     i  • 

S*  1   '.;    '     1 

DITOP*  R6.3  '  )  ) 

TS7  H»  r  ','>■■■',  .    '  ,•• 

T ST S«  ]  -  1  ?  /  .  •'     1 

STCRIT"  '--•.....,» 

TS"H"  Rsy^n. 

TSBS"  ?29  )]  .   i 


END 


;r    VIO     nPF  -,;  ■  t  T  '  0     .vsTt"     y/F^j'"       1       'rVl«5TPN 


l/«H/7<»     GEMt'RATI 


V3  '/"M 


JOB 

HAn  Dl  T           rHA     ••,[■ 

.  ^=5 

/     JC 

.  V- 

13^ 

L I \ r  c    •■'■!.'       -i 

-      1  i?e 

/     K     L»i 

•  1  f 

88 

C  A  p  n  F     •  -  "  • 

a         -1 

/     <     C  n 

•  1  ": 

ef 

P  L  0  T  T  ;   -     ,  •  -  r  i  r  n  s 

S         .  ~c 

/     1 

•      ! 

15 

HOnEI             SKCC.'JPS 

rv;"-  .  - 

/     MP     -• 

•  i 

PC 

HODEI      >"      -    CCNT  S 

5  1  c   •  r 

/    mO,,  ,.-. 

•  1 

TCT*L 

rH     -     ._      .. 

•  /'- 

^EIR    <*9?     1--'J1        LFOGtp    CH-T     -  '   /  '  1  /  7  /,     22:  19. 


:■  6<  33    I..FT  r     AF  TEH? 


■    I     l_  n  i  I  !■  '- 


20^- 


As  we  see  the  results  obtained  for  the  500  ft  aluminum 
case  are  acceptable  and  agree  with  all  our  considerations  in 
the  derivation  of  the  formulas,  namely  T0P>B0T. 

The  following  two  runs,  for  500  ft  steel  and  1000  ft  alu- 
minum  show  that  the  results  are  in  disagreement  with  T0P>B0T 
used  in  the  formulas.  We  will  see  next  why  this  happens  and 
derive  the  proper  set  of  formulas  to  be  used  for  such  cases. 
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VIO  OPERATING  SYSTEM  VERSION  1  REVISION   C!l2  03/K*/7*  GENERATED   PI *♦  / 3 P) •  7 *♦ 
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Id  LINES  PRINTED  PR2  * 

£8  CARD?  READ  * 
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1  •  -t-AC*AE  ) 
OP 

10)  TqP 

1  '  ,'d:X,   ,TrP='FU'.4  ) 
2*  )6CT 
S  2  0  X  $  'BOT-'F10»<H 

he  spacing  ,-f  longitudinal  stiffening  members' 

3v*  )S 
1 ,  2?X, 'S«»F10«*) 

CP  +  6* ( TOP  +  p-T  ) 
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STREhGHT  OF  SHIFS-PN  Tint  NfUT'-AL  AXIS  LOCATION 
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VIC  OPERATING  SYSTEM  VERSION  l  REVISION   7)12  03/04/7^  GENERATED   0<t/3CV7<»   1 


ion 

3f 

Be 

?2 

'04? 


HANDLING  CHARGE  *   .35  /  jnB  .35 

LINES  PRINTED  PR2  •  1.25  /  K  |_N  .  16 

CARDS  READ  *  1.50  /  K  CO  *13 

PLOTTER  VECTORS  $  .?[•,    /    \     ?,  .%?, 

MODEL  73  SECONDS  S?5.00  /  HOUR  *15 

MODEL  ^  SECONDS  512.50  /  HOUR  .00 

TOTAL  CHARGE  *     *79 


EIR  490  14731   LOGGED  OUT  05/01/7*  88:26.   *   2^.8^  LEFT  AFTER 


03  LOGINS- 
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STEEL 

POI  =  0.3 

E  =  30*106 

AL  =  500 

B  =  AL/8  =62.5 

DRA  =  B/3  =  20.833 

D  =  AL/9=55.556 

HM  =  1.371^3*B*DRA*AL2  =  44642.187*10 

SM  =  B*DRA*AL2  =  32551 . 562*10^ 


4 


9  -6 

W  =  0.30396*(l-POI'i)*DSC/E  =  276.6*10 

z  =  60.129 

HTOP  =  8 
HBOT  =  20.833 

STR3T  =  0.152222*HTOP*Z  =  4402.865   note  that  this  value 

should  be  corrected  for  the  Z 
approximation  and  actually  is 
slightly  smaller 
STR3B  =  0.152222*HBOT*Z2  =  11465.61 
We  see  from  the  tertiary  stresses  obtained  that  we  end  up  with 
the  following  diagram  of  stress  distributions: 

20597  -^0:L 


/ 


/ 


11466  U — 18534 
30000     ' 


/ 


/ 


/ 


N.A. 


25000 
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From  the  diagram  we  see  that  theneutral  axis  will  end  up  closer 
to  the  bottom  and  so  we  must  consider  BOT >  TOP  in  the  formulas. 


The  same  type  of  diagram  will  be  obtained  for  1000' aluminum  case 
since  the  tertiary  stresses  obtained  in  such  case  are 

STR3T  =  1498.73 

STR3B  =  9873.078 

For  the  1000'  steel  case  the  same  type  of  diagram  will  occur 
due  to  the  tertiary  stresses  being 

STH3T   =   4402.86 

STR3B  =  29004. 
Besides  this  situation  requiring  us  also  to  consider  the 
neutral  axis  closer  to  the  bottom. still  another  problem  arises 
since  STR3B  is  larger  than  DST=25000.  This  case  is  solved  by 
properly  correcting  the  inaccuracy  that  was  introduced  by  the 
Z  approximation  as  exemplified  next. 
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Now  let's  derive  the  formulas  capable  of  handling  the  500' 
steel  case,  the  10001  steel  case  and  the  1000'  aluminum  case. 
We  must  consider  BOT  >  TOP  and  DIT0P>DIB0T 

DITOP  =  D*(B+D)*B0T/P 

DIB0T  =  D*(B*T0P+D*B0T)/P 

I  =  ^*D2*B0T*( (3*B2+2*B*D)*T0P+(D2+2*B*D)*B0T)/P 
where 

P  =  2*D*B0T+B*(T0P+B0T) 
Starting  again  with  the  equation 

STH*DIB0T  =  SBH*DIT0P 
and  making 

A  =  DST-STR3T 

AB  =  DSC-STR3B 
we  obtain 

A*DIB0T  =  AB*DIT0P 
Knowing  that  now  the  top  is  the  thinner  plating  we  anticipate 

S  =  Z*T0P 
and  may  write 

STR3T  =  0.152222*HTOP*Z2 

STR3B  =  0.152222*HBOT*Z2 

S 

where  the  Z  value  in  STR3B  should  be  which  has  a  smaller 

BOT 

value  than  Z. 
We  then  obtain: 

A*D*(B*T0P+D*B0T)/P  =  AB*D*( B+D) *B0T/P 


21** 


Simplifying  and  making 

AC  =  (AB*(B+D)-A*D)/(A*B) 
we  obtain 

TOP  =  AC*BOT 
Flaking  now  the  top  to  be  stressed  the  most  in  tension  we  have, 

STH  =  HM*DIT0P/I  =A 
Replacing  DITOP  and  I  by  their  equations  we  have 

HM*D*(B+D)*BOT 

A  = 

*J.*D2*B0T*  (  (  3*Bb+2*B*D )  *TOP+  ( D2+2*B*D )  *BOT  ) 

and  obtain 

HM*(B+D)        3*B2+2*B*D 

BOT  = *TOP 

A*4*D*(D2+2*B*D)     D2+2*B*D 

Making 

AD  =  HM*(B+D)/(A*4*D*(D2+2*B*D) ) 

AE  =  (3*B2+2*B*D)/(D2+2*B*D) 
the  formulas  become 

BOT  =  AD-AE*TOP 
and  substituting 

TOP  =  AC*BOT 
we  finally  have 
AD 


BOT  =- 


1+AC*AE 
The  results  obtained  for  the  examples  chosen  are  presented  next: 
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We  could  have  obtained  TSBH  =  30000  by  correcting  the 
approximation  used  Z  =  S/BOT  as  is  done  for  the  1000'  steel 
case  that  follows  next. 

In  this  case  we  already  anticipate  the  use  of  a  Zl 
value  as  done  in  the   previous  section  and  according  to  the 
output  obtained  we  will  then  obtain  the  correct  result. 
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Now  we  see  that  with  an  assumed  value  Zl  =  35  we  obtain 

S  =  215.1808 

BOT  =  3.686        so  Zl  actual  =  58.38 

TSBH  =  47512. 793  (too  high) 

TSBS         also  too  high 

The  Zl  actual  value  obtained  suggests  that  we  should  have 
assumed  an  higher  Zl  value  and  with  an  assumed  value  Zl  =  40 
we  then  obtain 

S  =  203.626 

BOT  =  4.428        so  Zl  actual  =  45.98 

TSBH  =  34131.855  (too  high) 

TSBS         also  too  high 
This  suggests  an  assumed  Zl  value  higher  than  40  but  now  we  may 
interpolate  with  the  two  results  already  obtained  and  determine 
the  Zl  value  that  will  bring  TSBH  =  DSC  =  30000  psi  and  done 
properly  this  new  Zl  assumed  will  end  up  being  Zl  actual. 
From  interpolation  we  obtain 

(40-35)*( 34131. 855-30000) 

Zl  =  40+ =  41.544 

(47512.793-34131.855) 

and  now  assume  this  value  obtaining  the   following  result: 


226 


STReMGhT  Of   '•  I      '   'I  '  r    EUTRAj   AXIS  LOCATIPN 


f  AGE 


r 

r 


•  A  L  *      I  '-     r  H  F      Ll!"   T 
,"  L  =  1  .  •  •     . 

i  =  t  I  /  '• .  t  <-- 

p  =  A  L  / 1  *  • 

•  p'  1^-     I'       rHI      nuISSGM'q     <ATIO    HE     IMF    MATERIAL 
E  U  K    e;  T  E  K I 

p  0  I  =  (     .     J 

*u*     T  r-      '     E     '<  '  ,  K     iiS    virU.ijS    OF     THE     'ATERlAL 
i.  01 >    F  t  [  • 

E  *  3  (   •  *  !     i  *  *  i 

»051  x      I        TE'l         i-'i       '  Tits'?     I"!    tEKStQN 

h  S  T  =  ? ! 

*t'SC"      '        THE     nr«J!'jM    StRESS    IN    COM  Pi- ESS  ION 

•  --  r.  =  >i.  ■ 

>  l  ?  i     I '      M.i      >■'{;'■■■    :-■     ■  i  r       pi-     if  "tf-  T 
li'  =  1  .  3  /  i  .,  3  «  I    »M.'      -ai    *  •  | 
■  ■-  I    .      !  '■:      r  I  -       ■  ,  r;     '    ' ".        ;  I  i    t  \ \t     f^of^Eh  T 
>.;-  =  ■    •  ■■■>-      ■ ;  |   *A| 

!  =  l  .  :.-■    3:  (.',*  (  1  •  -E"".'  f  *cn  I  )  *L'sc/E 
i,=sprt<     ) 
•=i  ./U 

*HTCR*     IS    Tt,i-     h1-  '.  i.)    op    i,  -  r  c  f-'    Dm    DE  Cv 
Ml  of  =-  . 

*>RCT*     -T '-'     TEfE       EiH'Ap     if-     k'ATER    3V    BOTTOM 
MBCT=OR 

•  1=4  ■ .  4 r  .  */,  i  ?i  ,  rr->?/i  »-3£.    .  -''•  - 

S  n    P3       =1.1    5  ?  ?  ?  1  *  '     '•  ''  T  *  .'  T   *  Z  1 

A  =  )ST -S  T  E 3 T 
AP  =  DSC.  '-  Tf •  ~-t 
kCs(  4R»  (  p+p  )-a*i>  )/(  a  *n  I 


>  - 


+  0  )  /  ( 


*L,*(r,»U+2,+B*p)J 


.'\ r-  -  (   1  •  >-     »  ■  +  c  •  *  n  » ,"■  )  /  (     „  i  -f  ,_> ,  »  p  *  p  ) 

KTCP*     t'      THE     TMlfKi.'ESq    '  E     THE     TOE'     PLATING 

*\  f.  l  *     IS     t  k  E     THlC<^E3t!    OF    THE     BOTTP^    Pl.ATUn 

<  '.!  =  /.  ■  /  t  i  .*tr,^   , 

rOP  =  aO"OT 
iPlTf  IbMi   I  Top 
If    EC** '"•ATI  '  1 '  #2Zy,  'T3P='PJ     .  '♦  ) 
«RlTf  (  ••.  "V  )i  --t 

♦  S.     IS     I  EiE     SPArI     G    nr     L'    jGITUOIMAL    STIFFENING    KEMBERS 

S  =  Z  *  T  ")  t> 

.-i-'ITt-  (-,-...   |S 
3f      pOF^flT  (  '      •  ,  ??Y  ,'-;=•  r  j  .-  .  4  j 

p  =  p.»ii.,   r;t;   *,jrc,+QnT] 

#01  TOP*   Tn  THC  rr^^A  re     "IF  TMr  T[^p  rcnv'  THE  ^EITRAL  AX'S 
0 1 T  OP  =  T   (P*BOT  +  r*Wp-T)/P 
WITE(K/*  I  01  tv 

U(   '•  '  F  I  ..  T  (  '   •  ,  i  r  v  ,  •  -  I  t  o  - ■  -  «  E  l  •  ■ .  M 


227 


ST  RF  Nfi  h  7     0 >:'    S  P  1  f- c-     n  N 


C  ~  ^  '* 

i! :-  >•  <♦ 
>?>  ?-  <j  r 

I'.  ™  ,* .  ■ 

I     ! 

It  If 

I     !  ! 

I  ti»  i 

I  fc  J  ,  , 

[     1 

1*24 

fe  '-  .;  - 
1*7* 
It  Ql» 

I  •.  A  E 

I  *  C  A 

1 1 . '-  4 
ME      • 

l<  [I 

|P5  • 

I r- 1 1  << 

*.:   !    u 
tti '''  p  ' 

■  '        <:- 

>  t| 

wee 

0'   7- 

I    [  s  ]     .|i 

jf CV]      Ppj 

'-  CV] 
CV] 
ttC  CV  ] 
.    CV] 
CL1 
RCLJ 
CV] 
[V] 
W  CL] 
^  V  r. 

F 


^' 


P,f 


1   "' 

IK 


7 

s  t  r  3  a 

(  F 
2(? 

T^TP 

TsrH 

lia 


,I)T  r  ,-t  , 
■  I  I   <"T  =  ,    • 
.•K  J  T F  ( e  . 
-  '  ,  »   \T  f  i 


rhf 

1  hr 


JTRA|       AxIfJLnCATIf.Ki  p» 

i)F     Tur    BOTTOM    FROM     THF     NEUTRAL     Ay  Is 


r,  '■  I  =  i, .  ♦ ! 

*  TSTH*  IF 
*TSTS*JS 
*TSpr  »  I  S 

.  T  c  (  s  *  I ! 

*  S  1  f  :   r T  - 

*  a  ?■•  c  R I  t  * 


■  » i 

•r  ♦  ( 
i .  I 

1  Hf 
ThF 
TiJf 

T  C 

T  C 


T  S  i  '  ■  =  >  * 

t  a  t  s  s  a  ■ 
sT i  r i r  = 

I  ''.n  t  ■- 
f  1 1 1  ■     r  i 

M  K  I  7  E  (  p 
'  M  A  T  | 

'TT£(  - 

r OPf   •  T  (  ' 

X  S  ■ : '  ■  =  H  ■  •  < 
T  S  •  $  =  c  i  - 
■  C  '■  !  T '  =  I 
IkK  I  T  F  ( '  p 
F  (  f-' '  a  T  (  ' 
•  ^  I  T  p  (  -  , 

pOP?r/\  t  (  '     > ,  i  <•..,. 

rt^ITE  (  ^  ,  i  1  /  )  c., 
FORf  .'.  T  ,  •      •  ,  i  7v 

p-Nfl 

.-  A  7  C  C  v  ]  i 

^6 A  4  [V]  E 

P6CCCV]  <5M 

«6£i.  [V]  HTOP 

V  7  H    [  v  J  A 

ylhi   Cv  ]  PC  1 
.;•  75(   [  V  ] 

P»75C  CV]  DlROT 

,,,  7  7^-  [V]  STcRTt 

077C  r  V;  TSfc>S 

l  l '  I r-  i  .v 


TS     T'-r     '    1ST  a 

(  ;-  *  TOe  +     •  ><r  t  )  ,i 

•--•)[    i  .  ,    - 

1  «  v  ,  '     rFnT-'i-ipi.it) 

k/p^-,7     Or     I\F'-'TIa    of     tHf    \    SFrrinr- 
- . - ■   *  r  +;   ,*D*D)*70P*«OT+(0*D  +  2.*B*r 

^t1  '  ^c    '  •     i"qp  nuc    rn   i-  no^iNJG 

<"!'-  '■    :c     '         TOP    DljF     Tf      LAGGING 
^T'-:';r     t        Bqt    ?L  F     Tr^    (-OGGI^G 
PT^i  q?     :  ■      BOT    01  ■?     TO    <    \GGING 
Tur-    C^lTlpAL    STRFSS    rN     TOP 
Tur    CRITICAL     STRESS    '^     BOT     • 
►OlTOP/q      T  +'•!'-  tT 
1    ITUP/';  *  J  +f  T^oT 
1    iR*TOp*>/(  f  ,3     3 

■  *>  /  )  T  S  y  M 
'      •  /  l^v  j  •  TSTM=iF  lvl.3  ) 

•  '        )   T  ^  T  Q 

1      *  j  \Hy,  i  TSTc=,Fl  1,1.3  ) 

>■-•)  si  »p  j  r 

■  • ,  i.7v,  •s-TrpjT='Ftr.3) 

" '.  1 1  CT/i:     t+r'.lc;8?P2*WBn7'»F*c,/( 

•  '  ■  1  ■  '-  T  / ,  ■     l  + 

IT*!    Hi    .  -  /   ( 

'  ,  1«^,  •  TSF  M=.F'1C!.3) 


A  ^0  it 


•  A  . 
T  ,  /P 


oi *pni  )  > 


.2^32  *WE*0  r 


•■  *  S  /  I 


F)i    T*Bf  T  ) 
B  0  T  *  R  r  T  ) 


'3?6*?*S*(l»-r0I*POI)) 


,  'TSF 


=  i  (•  1  ;• .  3  ) 


1  95c*G  Its  '  Fl 


3  ) 


•<>>/.  Cv]  R 

1  '■'  ('  0  C  c  ]  .  P 

■6pprv]  •„ 

tfrC\/]  HBCT 

■71  /,  Cw]  SB 

7/fCCv]  TOP 

2AC  C|  ]  "5? 

3 ."  v  '  i  ]  5  '7 

1+  M  /,    '  i     ]  A  ? 

7,:'  [v]  SBf  RlT 


(  6 !  C  C  V  ] 
c  a  r  '♦  r  ;  ] 
r  a  r  ■  C  C  V  ] 
P  h  r  "•  C  V  ? 
t  7i  X  CV] 
(•  ?  '  '»  r  L  ] 

t    7  cr  *♦    C  V  ] 

f  7^  •'  CV] 

t'v  a  CL] 
F5f  pi  rL] 


DP  A 

?!  ft  9  <♦  C  V  - 

r 

PS7 

pi  (y  R  c  C  v  • 

r^    c 

U 

■^OOC  [S- 

c  --  ■_-  r 

S  T  F  3  T 

"i  ft  F ;   c  V  - 

71 

AC 

r.-  7  p  f;  C  V 

/I    - 

10 

pejac--  ca- 

&  I 

P 

P758 CV* 

nn    P 

Gf-'I 

^  7  6  C  C  V  • 

T^IT- 

72 

R50P  [L- 

R    ' 

9C 

P1626  CL- 

] 

Ir^PAM    |_APpL 
2  V)  7  "21     *  v  t,  I  • : » 
-•  E  .'  2     .  7  E  P  0 
EC6A     A I NJT 
rD04     %a 


2FC/!  •>. 

2  0  ft  (  t-t 

2  8  F  t  ~ '   f '  T 

2  f  r  (  .  i 


PC 3A  t' 

'-  L"  A  2  •  R  A  «  Q 
2  E'OC  •  1 1 


290C  A  LOG 
2F'5F  ,MES 
20PF  .5 


4 


3  4 


2P1P  «I 
2CPC  • W 

PATi  a EXP 


P7^  4 
2  A  E  A 

2E«  t 


S"RY-P0iNTg: 


228 


Tl      * 

3.3     35 

hO  7  = 

l».7**!    <* 

"l  ~ 

1    <-l  V         ,       >0 

n  i  t  o  => = 

hf    .  ,■•  .3  1    J 

DTH      T  a 

31  .3'''    ^ 

-  --. )  w  = 

;.>. c,t  ■) . 

Tq« 

i  p^a  1 .  - 1  -s 

'-  i  i  !•  :  r  = 

?999<s.o.   ? 

i  •--     i = 

3',     1  ?,  1   .  <     !|   ? 

T'             - 

?u^ u  .•+--<* 

'  >,;■•■  IT" 

6 1  9  *  -> ,  <  ?  3 

C!:     V  10     ^oFuirf.iG     '-  .  r  -, 


\/EF  c|nfj     i     REVISION       '•  l; 


3/04/74    uE UPRATED 


i  /  ?  ?  /  7  4 


F!  HA>  r)L  I  N'n    r'."' v    ■' 

:■  1  L  I  N  ES    P o  I N T  F  i      •       | 

:-  °f  r  A  ~  n  S     '•  f  A  ! 

I  '"  dL'-tTf  ^    wLir-s 

?"  -■• : i ^  c  i     7 , -.    s  t  r  r 

[  •/  y  0  n  r  [      r  ci     c  E  C  C  '* 


•  ~->r> 

/ 

J  r  ; 

J      3  .  ?? 

/ 

K     [_ . 

1       1   .  S  : 

/ 

i<       "  -^ 

•  -i  r 

/ 

1  P  ■  ■-> 

"-  /•   &  •     •  ', 

/ 

l  n   in 

512.  '-■:.-- 

/ 

i    n      -. 

TD1  u 

*' 

|  t  K      .- 

•  3C 

•  1  A 

•  13 

•  &<■> 

.   1   K 


IP     49P     1^731        L  it  H    C'UI     7-?/'o/7i 


^3; PI  • 


•  87    L  L  r T     AFTER 


fe    L     i  I l>  S 


229 


We  see  that  TSBH  is  just  a  little  too  high  and  the  same 
happens  with  TSBS;  here  it  might  be  the  case  that  TSBS  should 
have  been  the  one  under  consideration  for  the  choice  of  a  Zl 
value.  Nevertheless  for  the  purpose  of  this  work  this  example 
shows  the  procedure  to  be  used  along  with  the  equation  derived. 
The  accuracy  obtained  can  be  total.  Here  we  assumed 
Zl  =  *U.5^4  and  ended  up  with  Zl  actual  =  41.83  which  is  close 
enough  for  our  purpose  and  TSBH  and  TSBS  is  very  close  to  what 
it  should  be  to  use  the  material  effectively. 
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This  last  example  shows  that  TSBH  is  low  and  the  Z 
approximation  used  in  the  formulas  here  again  should  be 
corrected  as  already  exemplified  in  order  to  bring  it  closer 
to  the  DSC  value. 
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CONCLUSION 

From  this  work  the  author  is  led  to  conclude  that  there 
is  no  simplified  approach  to  the  neutral  axis  location  in  the 
midship  section,  capable  of  handling  properly  all  or  only 
a  majority  of  cases.  It  is  advisable  to  look  into  each  parti- 
cular situation  by  first  gaining  some  insight  at  what  the 
diagrams  of  stress  distributions  may  look  like  as  exemplified 
in  several  sections.  The  two  last  sections  exemplify  well  why 
this  may  be  of  some  advantage  and  specially  the  last  section  which 
shows  that  while  at  first  it  was  expected  for  the  examples  stu- 
died a  neutral  axis  closer  to  the  top,  actually  the  results 
ended  up  requiring  a  neutral  axis  closer  to  the  bottom  and  a 
new  set  of  formulas  had  to  be  derived. 

As  a  future  development  for  this  work  it  is  here  sugges- 
ted that  a  more  ship  like  cross  section  could  be  analysed  as 
a  step  forward  towards  obtaining  a  method  of  determination  of 
scantlings  as  direct  as  possible.  This  next  step  should  consi- 
der the  derivation  of  proper  formulas  for  DITOP,  DIBOT  and  I 
involving  only  the  plate  thicknesses  as  unknowns.  It  should 
be  considered  then,  in  order  to  simplify  the  problem. ,  a  stan- 
dard type  of  frame  and  its  scantlings  should  be  constant  throu- 
ghout the  midship  section. 
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